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SUMMARY 


This  work  contributes  new  understanding  of  the  dominant 
processes  for  phase  transformations  which  may  occur  under 
conditions  of  shock  wave  loading  of  materials.  It  should  be  of 
interest  to  material  scientists  and  design  engineers  who  deal 
with  shockwaves  or  impulsive  loads.  The  work  was  supported 
principally  by  the  Naval  Surface  Weapons  Center  Independent 
Research  Program,  Task  177,  with  additional  support  from  the 
Air  Force  Office  of  Scientific  Research  under  Grant  No . AF-AFSOR- 

1088-66  and  (ARPA)  Contract  F-44620-67-C-0087.  - 
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CllAPTKR  1 
INTRODUCTION 

UnOoi  tlio  InlliUMico  ol  .shock  waves  certain  materials 
iindei>.jo  polymorpluc  t lanstoi  mat  ions . The  time  available  in  a 
:'.lK)ck  experiment  tot  siicli  a t rane.tormat  ion  t>.>  occur  does  not 
exceed  a lew  microsecond:;,  yt't  t lu'  ::ame  t ranst\)rmat  ion  may  re- 
quire minutes  oi  houi  lor  completion  in  a :;tattc  experiment. 
This  larqe  t ran:;l  ormat  ion  rale  i:;  an  important  di:;tinctive  tea- 
ture  ot  I ran:;l  ormat  ions  in  sluH'k  wave.s.  Little  is  known  about 
atomic  mechanism:;  which  cau:;e  these  ultra(a:U  t ransl  ormat  ions . 
Thi.s  void  in  knowledqe  is  partly  ilue  to  lack  ot  experimental 
data.  The  ne<Hl  lor  data  and  uiuler :}t and i nq  ot  the  t ranstormat  ion 
proce:;:;  is  the  rea:;on  toi  I hi:;  work. 


l.l.  General  hack'jro;ind 
of  Polymorphism 


The  ability  ot  a compound  to  ctystallize  in  ditterent 
structures  is  we  I l-e:;t  abl  i shed . hri;himan^  showed  that  many 
materials  underqo  a t ; anst ormat ion  to  a new  structure  under  tiu' 
Intluonce  of  pres:;ure.  In  more  rec«'nt  times  it  ha:;  been  shown 
that  many  tran:;tormat  ions  ob:;erved  :;l,;tically  ai  e also  observed 
dynamically;  i.e.,  under  shock  comt're:;;;  ion . Kxamples  ot  this 
are  found  in  materials  such  as  iron,''*'  carboi;,  bismutl;,^'^ 

sillcon,^'^  boiani  nitride,  m.my  otlieis. 
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To  understand  the  varied  aspects  of  polymorphism  it  is 
necessary  to  relate  static  results  to  dynamic  measurements. 
Caution  must  be  exercised  in  making  this  relation  because  of  the 
distinct  differences  in  compression  methods.  A plane  shock  wave 
results  in  uniaxial  compression  at  very  large  strain  rates  accom- 
panied with  large  shear  forces.  These  properties  can  result  in 
defect  formation,  plastic  flow,  and  other  irreversible  effects. 

Some  progress  in  relating  static  and  dynamic  results  has 
been  made  by  identifying  shear  deformation  as  a likely  mechanism 
that  accelerates  transformations.  The  addition  of  shear  to  a 

sample  under  static  compression  has  been  shown  to  accelerate  the 

12 

transformation  process.  The  accelerating  effect  of  adding 

shear  to  static  experiments  has  led  Leiserowitz,  et  to 

draw  a controversial  analogy  between  experiments  on  shock  wave 

compression  and  experiments  of  hydrostatic  pressure  with  shear 

stresses  added.  The  role  shear  deformation  plays  in  accelerating 

transformation  rates  in  dynamic  experiments  is  difficult  to 

ascertain  since  an  elevation  of  hydrostatic  pressure  above  that 

normally  required  to  produce  the  transformation  can  greatly  in- 

14 

crease  the  transformation  rate.  To  date  no  explanation  for 

differences  between  static  and  dynamic  rates  is  satisfactory. 

In  dynamic  experiments  the  transformation  rate  affects 

the  shape  and  amplitude  of  shock  waves  which  evolve  towards  the 

steady  state.  Studying  evolution  of  a wave  in  time  provides  some 

understanding  of  dynamic  polymorphism.  Experimental  shock 

studies  of  the  kinetics  of  solid-solid  phase  transformations 

2 15 

have  been  limited  primarily  to  four  materials:  iron,  ' 
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16  17  18  19 

potassium  chloride,  ' antimony,  and  Arkansas  novaculite. 

A framework  for  interpreting  experiments  in  dynamic  polymorphism 

20 

was  provided  by  Hotie  and  Duvall,  who  combined  continuum  me- 
chanics and  thermodynamics  to  treat  the  evolution  of  shock  waves 
in  materials  with  time-dependent  phase  transformations. 

1.2.  Alpha  to  Epsilon 
Transformation  in  Iron 

Discovery  of  and  serious  attempts  to  understand  the  iron 

21 

transformation  appear  to  have  followed  the  publication  by  Walsh 
of  high  pressure  iron  Hugoniot  data  which  disagreed  markedly  with 
existing  work  on  its  compressibility  at  lower  pressures.  In  par- 
ticular, the  results  implied  compression  of  iron  significantly 
greater  than  would  be  expected  from  extrapolation  of  static  com- 
pressibility data  of  Bridgman.^  This  stimulated  studies  using 
shock  waves  at  lower  stresses.  The  discrepancy  in  compress- 
ibility was  resolved  by  the  discovery  by  Minshall  and  his  co- 
2 22 

workers  ' of  multiple  shock  stress  waves  from  which  they 
inferred  the  existence  of  a polymorphic  transformation  near 
130  kbar. 

Many  questions  still  remain  unresolved  about  the  pressure- 
induced  iron  phase  transformation.  The  mechanism  of  transforma- 
tion and  effects  of  shear,  temperature,  pressure  gradients,  and 
strain  rate  on  the  transformation  process  are  not  fully  under- 
stood. In  addition,  Hugoniot  data  and  equilibrium  thermodynamic 
predictions  differ  for  the  mixed-phase  region.  The  work  in  this 
dissertation  is  directed  toward  understanding  the  non-equilibrium 
behavior  of  the  mixed-phase  region  for  the  iron  transformation. 
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1.2.1.  Background 

Iron  is  an  elastic-plastic  material  which  undergoes  a 
phase  transformation  when  stressed  above  130  kbar.  These  prop- 
erties will  result  in  three  forward-facing  shocks  (elastic, 
plastic  I,  and  plastic  II,  respectively)  traveling  away  from  the 
impact  surface.  Under  equilibrium  conditions  the  stress  behind 
the  plastic  I shock  is  a measure  of  transformation  pressure  for 
the  polymorphic  transformation.  Time-dependence  of  the  transfor- 
mation requires  that  stress  in  the  plastic  I shock  and  the 
plastic  II  shock  overshoot  their  equilibrium  values  for  short 
times  after  impact.  The  evolution  of  the  shocks  is  illustrated 
in  Fig.  1.1  where  shock  wave  pressure-distance  profiles  are  drawn 
for  two  different  times  following  impact  on  the  left.  Measure- 
ment of  shock  amplitudes  as  a function  of  thickness  can  be  used 
in  conjunction  with  appropriate  models  to  obtain  information 
about  kinetics  of  the  transformation  process. 

1.2.2.  Shock  Profile  Data  Related 
to  Transformation  Kinetics 

Transformation  rate  affects  the  shapes  and  amplitudes  of 
shocks  in  iron.  There  are  four  types  of  existing  experimental 
data  on  shock  profiles  in  iron  which  relate  to  transition  kinet- 
ics. The  four  types  of  data  are;  (1)  amplitude  of  stress  behind 
the  plastic  I shock  as  a function  of  sample  thickness,  (2)  ampli- 
tude of  stress  behind  the  plastic  I shock  as  a function  of  final 
driving  stress  for  constant  sample  thickness,  (3)  rise  time  in 
the  plastic  II  shock  front,  and  (4)  thickness  of  the  plastic  II 
shock  front  inferred  from  residual  metallurgical  effects. 
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FIGURE  1.1  TWO  SHOCK  WAVE  STRESS-DISTANCE  PROFILES  IN  IRON  FOR 
TIMES  tj  AND  WHERE  <t3. 
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The  plastic  I wave  stress  exceeds  its  eventual  time- 

independent  value  near  the  impact  boundary,  and  time  required  for 

the  stress  to  decay  to  its  time-independent  value  is  related  to 

2 2 

kinetics  of  the  transformation.  Minshall  and  his  co-workers,  ' 
measured  a slow  decay  in  the  plastic  I wave  stress  as  a function 
of  thickness  for  a decaying  driving  stress.  The  plastic  I stress 
was  about  140  kbar  in  a sample  6-mm  thick;  it  decayed  to 
130  kbar  in  a sample  40-mm  thick.  In  the  present  experimental 
work  I measured  a similar  slow  decay  in  plastic  I stress  as  a 
function  of  thickness  for  constant  final  driving  stress.  The 
plastic  I stress  was  about  139  kbar  for  a 1-mm-thick  sample 
and  decayed  to  131  kbar  in  a 25-mm-thick  sample.  Final 

15 

driving  stress  was  near  200  kbar.  Barker  and  Hollenbach 
measured  a small  increase  in  the  plastic  I stress  with  increase 
of  final  driving  stress  for  constant  sample  thickness;  amplitude 
of  stress  behind  the  plastic  I shock  increased  by  6 kbar  when 
final  driving  stress  changed  from  130  kbar  to  300  kbar  for 
6 . 35-mm-thick  samples. 

The  transformation  process  occurs  principally  in  the 

plastic  II  shock  front;  therefore,  the  rate  of  transformation 

may  contribute  significantly  to  rise  time  of  this  wave.  I and  a 

15  23 

number  of  experimenters  ' have  measured  rise  times  of  0.2 
to  0.3  Msec  for  the  plastic  II  shock  front  in  iron  by  monitor- 
ing surface  motion  of  shocked  samples.  A less  accurate  method 

of  obtaining  rise  time  information  is  from  residual  metallurgical 
24  25 

effects,  ' which  imply  shock  front  thickness;  for  an  approxi- 
mately steady  wave,  rise  time  is  obtained  by  dividing  shock  front 
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24 

thickness  by  wave  velocity.  Smith  has  suggested  that  the 
plastic  II  shock  front  width  in  iron  is  less  than  0.02  mm, 
based  on  metallurgical  data.  This  shock  width  implies  a rise 
time  of  about  0.005  psec. 

To  interpret  some  of  the  above  experimental  results, 

20 

Horie  and  Duvall  developed  a model  of  mixed  phases  which  con- 
tained a constant  relaxation  time  parameter.  This  relaxation 
time  has  to  be  evaluated  from  experimental  data.  In  considering 

available  data,  they  found  a relaxation  time  of  1/3  psec  based 

23 

on  rise  time  measurements  of  Novikov,  et  for  the  plastic  II 

shock  front  of  0.01  psec  based  on  the  plastic  II  shock  front 

24 

width  of  less  than  0.02  mm  reported  by  Smith  , and  of  20psec 

based  on  the  slow  decay  of  the  plastic  I stress  observed  by 
22 

Minshall.  A relaxation  time  less  than  0.05  psec  is  required 

to  explain  my  data  on  decay  of  the  plastic  1 stress,  while  Barker 

and  Hollenbach^^  found  that  a relaxation  time  of  0.16  to 

0.18  psec  was  required  to  explain  their  data  on  plastic  I stress 

decay.  The  disagreement  among  values  of  relaxation  time  impugns 

20 

the  validity  of  the  Horie  and  Duvall  assumption  of  constant 
relaxation  time. 

1.2.3.  Evidence  for  Non- 
equilibrium Behavior 
of  Mixed  Phases 

Experimental  and  theoretical  evidence  exists  which  sug- 
gests that  in  the  pressure  range  of  130  to  300  kbar  iron 
is  not  in  equilibrium  and/or  the  transformation  does  not  go  to 
completion.  The  more  pertinent  results  are  summarized  in  this 
section. 
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At  least  two  theoretical  attempts  have  been  made  to 

understand  the  iron  Hugoniot  data  between  130  and  300  kbar. 

26 

Duvall  and  Horie  demonstrated  that  the  calculated  slopes  of 
the  equilibrium  phase  line  in  pressure-temperature  space  dif- 
fers significantly  from  the  measured  slope.  More  recently, 

27  28  29 

Andrews  ' ' calculated  the  equilibrium  Hugoniot  for  iron 

and  found  it  to  differ  from  Hugoniot  data  between  130  and 
300  kbar.  The  conclusion  from  both  treatments  is  that  dynamic 
transformation  data  do  not  ag  ae  with  equilibrium  thermodynamics. 

Single  crystals  of  silicon  iron  have  been  subjected  to 
x-ray  diffraction  and  metallographic  analysis  after  exposure  to 
shock  waves  and  relief  waves  which  caused  the  formation  and  dis- 
integration of  the  epsilon  phase. As  a result  of  this  cycle 
of  transformation,  one  would  expect  single  crystals  to  be  trans- 
formed into  polycrystalline  aggregates.  The  diffraction  analysis 
shows  that  the  material  had  retained  in  large  measure  its  origi- 
nal orientation,  though  it  was  polycrystalline.  There  are  two 
possible  explanations  for  this.  The  first  is  that  the  alpha  to 
epsilon  transformation  and  its  reversal  are  crystallographically 
reversible  to  some  extent. The  second,  suggested  by  German, 
et  al. , is  that  the  transformation  does  not  go  to  completion. 
This  suggestion  appears  to  be  in  better  agreement  with  various 
experimental  results. 

Much  of  the  experimental  evidence  implying  non-equilibrium 
behavior  comes  from  the  results  of  static  isothermal  compression 
experiments.  Such  experiments  show  that  the  iron  transformation 
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begins  near  130  kbar  and  goes  to  completion  for  stresses 
greater  than  150  kbar . This  behavior  disagrees  with 
equilibrium  thermodynamics,  which  requires  complete  transforma- 
tion to  occur  at  constant  stress  when  temperature  is  held 
constant.  This  requirement  appears  never  to  be  realized  experi- 
mentally. A very  thorough  investigation  of  the  Bismuth  I- 
Bismuth  II  transformation  by  Zeto  and  Vanfleet®  shows  that  the 
onset  of  transformation  in  their  apparatus  always  occurs  at 
greater  pressure  than  does  the  reverse  transformation.  They  and 
others  take  equilibrium  pressure  of  the  transformation  equal  to 
the  mean  of  direct  and  reverse  transformation  pressures. 

An  anomaly  exists  in  static  isothermal  compression  data 
which  has  not  been  carefully  discussed  in  the  literature;  In 
some  experiments  iron  appears  to  transform  completely  at  the 
transformation  pressure,  while  in  other  experiments  with  the 
same  equipment  the  lower  phase  persists  to  higher  pressures.  No 
attempt  is  made  to  explain  this  anomaly  because  it  involves 
analysis  of  static  compression  equipment  and  techniques  which  is 
beyond  the  ability  of  the  author  and  the  intended  scope  of  this 
dissertation. 

1.3.  Outline  of  This  Work 

The  basic  problem  studied  in  this  work  is  the  body- 
centered-cubic  (alpha)  to  hexagonal-close  - packed  (epsilon) 
transformation  in  iron  under  shock  compression.  Iron  samples 
were  all  initially  at  room  temperature  prior  to  shock  loading. 
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Measurements  were  made  of  the  evolution  of  shock  waves  with  dif- 
ferent sample  thicknesses  for  approximately  constant  final 
driving  stress.  Experimental  data  were  interpreted  with  the  aid 
of  a mixed  phase  continuum  model  with  constant  relaxation  time 
to  obtain  some  information  about  kinetics  of  the  process. 

Hugoniot  measurements  in  the  mixed  phase  region  are  shown  to 
correlate  well  with  a simple  relation  between  difference  of  Gibbs 
energies  and  fraction  of  material  transformed.  Some  factors 
relating  to  deviations  from  equilibrium  in  the  mixed  phase  region 
are  explored. 

1.4.  Summary 

The  reasons  for  choosing  to  study  iron  were;  (1)  there 

exists  a wealth  of  thermodynamic  data  on  iron,  (2)  effects  of  the 

time-dependent  phase  transformation  on  shock  profiles  have  been 

20 

theoretically  calculated  by  Horie  and  Duvall  and  later  by 
27  29 

Andrews,  ' and  (3)  the  elastic  precursor  stress  amplitude  does 

not  depend  on  final  driving  stress  for  sample  thicknesses  of 

, ..  36 

3 mm  or  greater. 

The  experimental  intent  of  this  study  was  to  measure  the 
evolution  of  the  plastic  I shock  in  polycrystalline  Armco  iron 
when  final  driving  stress  was  201  kbar.  Basic  findings  were: 

(1)  stress  behind  the  plastic  I shock  increases  from  131  to 
140  kbar  when  sample  thickness  decreases  from  25  to  1 mm, 

(2)  transformation  stress  measured  in  a 25.4-mm-thick  sample 
is  131.4  + 3.3  kbar,  (3)  little  or  no  variation  of  plastic  I 
wave  amplitude  for  propagation  distances  between  0.9  and 
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6.35  nun,  (4)  relative  volume  behind  the  plastic  II  wave  is 
Vj/Vq  a 0.871  + 0.008  for  a stress  of  201  + 8.4  kbar,  (5)  rise 
time  for  the  plastic  II  shock  is  0.18  + 0.02  psec. 

In  terms  of  the  Horie-Duvall  model,  decay  of  the  plas- 
tic I wave  implies  an  initial  transformation  rate  greater  than 
2 X 10^/sec  for  final  driving  stress  of  201  kbar. 
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CHAPTER  2 

EXPERIMENTAL  CONCEPTS 

In  this  chapter  a number  of  relations  are  given  for  con- 
verting experimental  data  into  more  useful  forms.  Sections  deal 
with  basic  relations,  wave  propagation,  and  free  surface  velocity 
measurements . 


2.1.  Basic  Relations 

The  equations  for  plane  one-dimensional  flow,  independ- 
ent of  material  properties,  are: 

H * • 12.1) 


if 

8x  ' 


(2.2) 


dE  _ 
■ 


and 


(2.3) 


where  t is  time,  x is  Eulerian  space  coordinate,  V is  spe- 
cific volume,  u is  mass  velocity,  E is  specific  internal 
energy,  and  P is  stress  in  the  x direction,  positive  in  com- 
pression. Thermal  conductivity  is  assumed  to  be  negligible. 

For  steady  waves  the  flow  equations  lead  to  the  following 
jump  conditions  for  a shock  discontinuity; 


j 

J 
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= Pj^{U2-Uj^)  , (2.4) 

P2  - Pi  = Pi(U2-Uj^)  (U2-U3^)  , (2.5) 

E2  - El  = j (P2+P1)  (Vj^-V2)  , (2.6) 

where  U2  is  shock  velocity.  Subscript  2 is  for  the  state  be- 
hind the  shock  front  while  subscript  1 is  for  the  state  ahead  of 
the  shock.  Density  at  room  temperature  and  atmospheric  pressure 
is  given  by  Pq.  The  measurement  of  any  two  parameters  of  the 
set  (P,V,E,U2»u)  for  known  conditions  ahead  of  the  shock  is 
sufficient  to  determine  the  other  three.  The  two  measured 
parameters  for  this  work  are  shock  velocity  and  change  in  free 
surface  velocity,  which  is  assumed  to  be  twice  the  change  in 
particle  velocity;  i.e., 

AUfg  s 2Au  . (2.7) 

37 

Walsh,  e^  al,  have  shown  that  particle  velocity  obtained  from 
Eq.  (2.7)  is  accurate  to  1 percent  for  a single  shock  in  iron 
with  final  stress  of  400  kbar. 

Free  surface  velocity  produced  by  reflection  of  a shock 
from  an  unconfined  surface  is  the  sum  of  particle  velocities  pro- 
duced by  the  shock,  u,  and  by  the  reflected  rarefaction,  u^ , 
respectively;  i.e.. 
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An  elastic-plastic  material  which  is  also  polymorphic, 
such  as  iron,  has  a complicated  release  path  in  the  pressure- 
particle  plane.  Rigidity  of  the  material  and  kinetics  of  the 
reverse  transformation  help  determine  the  value  of  u^. . Treat- 
ment of  the  release  of  stress  for  shocked  elastic-plastic  solids 
has  been  fairly  successful.  Measurements  of  the  reverse  trans- 
formation reported  by  Barker  and  Hollenbach^^  show  that  even  in 
this  case  Eq.  (2.7)  is  a good  approximation.  Equation  (2.7)  has 
been  used  throughout  the  present  work  to  convert  free  surface 
velocities  to  particle  velocities. 

2.2.  Wave  Propagation 

Plane  shock  waves  are  generated  in  the  iron  sample  by 
impact  of  a flier  plate  or  a plane  detonation  wave.  After 
traveling  a short  distance,  the  initial  shock  develops  into 
elastic,  plastic  I,  and  plastic  II  waves  described  in  the  intro- 
duction. These  proceed  through  the  sample  and  each  interacts  in 
turn  with  its  plane  free  surface.  A map  of  the  process  in  space 
time  is  shown  in  Fig.  2.1.  The  space  cordinate  is  Lagrangian, 
being  the  undisturbed  x coordinate  denoted  h.  The  Lagrangian 

I 

wave  propagation  velocity  is  denoted  by  a prime:  for  elastic 

If 

wave,  U2  for  plastic  I,  and  for  plastic  II.  Final  stress 
behind  the  plastic  II  shock  front  is  greater  than  transformation 

i 

stress.  Initial  time-dependent  effects  involved  in  the  formation 

! 

of  the  three  waves  are  not  shown  in  the  diagram.  Reflected  wave  i 

• - ^ 

velocities  are  denoted  R and  wave  paths  in  Fig.  2.1  are  ] 

I i 

labeled  according  to  their  slopes;  e.g.,  dt/dh  = l/U^^  for  the  'j 
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LAORANOE  DISTANCE 


FIGURE  2.1  LAGRANGIAN  DISTANCE-TIME  DIAGRAM  OF  SHOCK  WAVE  PROPAGATION  IN  IRON. 
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elastic  shock.  Elastic  reverberations  between  forward-facing 
shock  fronts  and  free  surface  are  ignored.  Times  labeled  1,  2, 
and  3 correspond  to  arrivals  at  the  free  surface  of  elastic, 
plastic  I,  and  plastic  II  shocks,  respectively. 

Point  4 is  the  intersection  of  the  forward-facing  plas- 
tic II  shock  with  the  backward-facing  relief  wave  originating  at 
point  2.  Point  4 lies  on  the  locus  of  h-t  points  which  nearly 
represents  the  path  the  original  forward-facing  plastic  II  shock 
would  have  taken  without  interference  from  the  elastic  relief 
wave.  Coordinates  of  point  4 can  be  calculated  from  the  results 
of  a series  of  experiments  with  the  same  impact  stress  but  dif- 
ferent sample  thicknesses,  h^.  The  cccr i t-.  . f voi.it  4 are 
given  by 


^4  ^ ^2 


(t3-t2)-U^ 


U,  + k 

c a 


(2.9) 


hjj  — hg  3 2 c 3 

I • 

u,.  + 

c a 


(2.10) 


where  hg  is  the  coordinate  of  the  free  surface.  Arrival  times 

I 

t2  and  t^  are  measured;  U^.  is  assumed  equal  to  elastic  wave 
speed  U|^. 

Eulerian  velocities  for  plastic  I and  plastic  II  waves 


are  obtained  from  Lagrangian  velocities  through  use  of  the  rela- 
.38 

. Pi 


tions" 


«2  = ^ ' 


(2.11) 


.! 
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U 


(U3-U2)  , 


(2.12) 


where  p and  u are  density  and  particle  velocity,  respectively, 
behind  the  wave  defined  by  the  subscript. 

I 

Elastic  wave  velocity,  = CJj^,  was  taken  to  be 
6.18  mm/psec.^°  Plastic  I velocity  (Euler ian)  for  individual 
experiments  was  calculated  by  assuming  the  wave  to  be  centered 
at  h=0,  t=0,  which  gives 


U. 


ho  + 


(2.13) 


Experimental  data  on  arrival  times  at  the  free  surface  for  dif- 
ferent hg  support  the  assumption  that  the  plastic  I wave  is 
centered  at  h=0,  t=0. 


2.3.  Free  Surface  Velocity 
Measurements 

Free  surface  and  shock  velocities  were  measured  using  the 
arrangement  shown  in  Fig.  2.2.  The  mirrors  shown  were  silvered 
on  their  lower  surfaces  and  collimated  light  was  allowed  to  fall 
on  the  assembly  at  near-normal  incidence.  Reflections  were  re- 
corded on  film  in  a rotating  mirror  camera.  Part  (b)  of  the 
figure  represents  the  film  record.  Portions  below  the  broken  and 
near-horizontal  lines  are  exposed;  upper  portions  are  unexposed. 

Shock  impact  against  a mirror  surface  causes  the  amount 
of  light  being  reflected  to  decrease.  The  change  in  reflected 
light  is  recorded  by  the  camera  as  a function  of  time.  Time 
events  labeled  on  the  record  are  for  flier  arrival  at  the  impact 


17 


FIGURE  2.2  reflected  LIGHT  PLATE  SLAP  EXPERIMENT  (a)  AND  CORRESPONDING 
STREAK  CAMERA  RECORD  (b). 
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surface,  plastic  I arrival  at  the  free  surface  of  the  iron 

iiample,  and  plastic  II  shock  arrival  at  the  free  surface, 

t^.  Elastic  shock  arrivals  at  the  free  surface  were  not  recorded 
because  of  their  small  amplitudes.  The  continuous  trace  due  to 
the  free  surface  impacting  the  tilted  mirror  has  two  linear  parts 
connected  by  a curved  part.  The  first  linear  part  represents 
free  surface  motion  produced  by  the  plastic  I shock  and  makes  an 
angle  with  the  distance  axis  on  the  film.  The  curved  part 
is  due  to  the  finite  rise  time  of  the  plastic  II  shock  front. 

The  second  linear  part  represents  free  surface  motion  produced 
by  the  plastic  II  shock.  It  makes  an  angle  Y2  with  the  dis- 
tance axis  on  the  film.  Tilt  of  the  flier  plate  is  indicated  by 
non-simultaneous  arrivals  at  the  outside  mirrors.  A,  labeled 

t 

tj^  and  tj^  on  the  record,  and  by  non-simultaneous  arrival  over 

I 

mirror  B on  the  sample.  Times  tj^  and  tj^  define  a line 
which  makes  an  angle  oj  with  the  distance  axis. 

Average  shock  wave  velocities  and  free  surface  velocities 
are  determined  from  measured  distances  and  slopes,  respectively, 
on  the  film  record  and  from  the  known  writing  speed  of  the  camera. 
A record  is  shown  in  Fig.  2.3  with  the  same  labels  as  in  Fig.  2.2. 

Film  records  of  the  free  surface  motion  were  read  with  a 
Vanguard  analyzer  or  a traveling  microscope.  Distances  on  the 
film  as  small  as  0.006  and  0.001  mm  can  be  resolved  with  the 
Vanguard  analyzer  and  microscope,  respectively.  To  accurately 
determine  slopes  of  linear  parts  of  the  trace, many  points  were 
read  and  fitted  by  least  squares  to  a straight  line.  Each  trace 
was  read  at  least  three  times  and  the  resulting  slopes  were 
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averaged.  In  order  to  preclude  points  of  the  curved  portion  from 

the  linear  fits,  the  entire  trace  was  graphed  and  limits  of  the 

linear  portion  were  determined  by  inspection.  This  judgment  was 

reviewed  and  proper  adjustments  were  made  after  least-squares 

fitting  of  the  straight  line  sections  was  completed.  The  time 

span  of  the  curved  portion  was  determined  from  the  record  and 

interpreted  as  a measure  of  rise  time  in  the  plastic  II  shock 

front.  The  equation  relating  free  surface  velocity, Uj2/  pro- 

. 39 

duced  by  the  plastic  I shock,  to  measured  parameters  is 


u tan  a 

_ c 

*f2  (tan  " tan  oi) 


(2.14) 


where  u^  is  camera  writing  speed,  is  camera  magnification 
(image/object),  a is  angle  between  tilted  mirror  and  the  sample 
surface,  y^^  is  the  angle  made  by  the  first  linear  trace  on  the 
film,  and  u is  tilt  of  the  shock  wave  measured  on  the  film. 
These  are  indicated  in  Fig.  2.2.  The  small  correction  due  to 
slight  rotation  of  the  surface  by  the  elastic  shock  has  been 
ignored  in  Eg.  (2.14). 

The  equation  relating  free  surface  velocity,  pro- 

39 

duced  by  the  plastic  II  shock  to  measured  parameters  is 


u^  tan 


f3  “ M,  tan  y. 


tan  CO 
tan  y2 


^ (^2~^f2\ 

yi  \U3-Uf3/ 

CO  /^2~^f2\ 

^2  V3““f2/ 


(2.15) 
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1 


where 


tan  a tan 

tan  ^ 


and  whci'e  Yo  is  defined  in  Fig.  2.2. 


(2.16) 


y 


I 

t 

I 


A 


I 
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CHAPTER  3 

EXPERIMENTAL  PROCEDURES 

3.1.  General  Facilities 

Most  experiments  were  done  at  the  Naval  Surface  Weapons 
Center  (NSWC) . The  facility  used  is  a test  chamber  with 
61-cm-thick,  steel-reinforced  concrete  walls  to  contain  the 
explosion  and  accompanying  noise.  It  was  equipped  with  a streak 
camera  and  a number  of  oscilloscopes.  The  Cordin  132  streak 
camera  is  capable  of  writing  on  the  film  at  speeds  up  to 
20  mm/psec.  Ancillary  equipment  for  calibration  and  for  syn- 
chronizing the  recording  equipment  was  also  available.  An 
explosive  preparation  facility  was  available  to  cast,  machine, 
and  press  explosive  charges  to  specifications.  A computer 
facility  consisting  of  a CDC  6400  computer  and  accessory  equip- 
ment was  available  for  use  in  data  reduction. 

The  preparation  of  target  components  and  their  assembly 
was  all  done  by  the  author.  Preparation  included  the  machining 
and  lapping  of  most  components.  Assembly,  preparation,  and 
handling  of  explosive  components  were  done  by  a technician,  as 
required  by  operating  procedures  of  the  laboratory.  This  techni- 
cian also  assisted  the  author  in  optically  aligning  test  assem- 
blies for  streak  camera  experiments.  All  other  aspects  of  the 
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experiments  were  done  by  the  author,  except  where  otherwise  indi- 
cated in  the  text. 

Two  quartz  gage  experiments  were  done  at  the  Shock 

Dynamics  Laboratory  at  Washington  State  University,  whose  facili- 

40 

ties  are  described  elsewhere. 


3.2.  Armco  Iron  Properties 

A 10 . 2-cm-diameter , 61.0-cm-long  bar  of  Armco  Magnetic 
Ingot  iron  was  obtained  for  this  program.  This  large  specimen 
was  obtained  to  insure  that  similar  test  samples  would  be  used 
throughout  the  course  of  the  study.  All  samples  for  this  pro- 
gram came  from  a 10-cm  length  of  this  bar  and  were  used  as 
received  except  for  cutting,  surface  grinding,  and  lapping  of 
sample  faces.  Metallurgical  properties  of  iron  and  procedures 
for  sample  preparations  are  given  below. 
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3.2.1.  Metallurgical  Properties 

Impurities  present  in  a single  unshocked  sample  of  Armco 
iron  were  obtained  by  the  NSWC  using  wet  chemical  and  spectro- 
chemical  analyses.  The  sample  tested  was  99.84  percent  pure 
iron,  and  accuracy  of  the  over-all  analysis  was  better  than 
0.05  percent.  Unmeasurable  traces  of  silicon,  nickel,  chromium, 
and  molybdenum  were  detected.  Measurable  amounts  of  carbon 
(0.14  percent),  phosphorous  (0.002  percent),  sulfur  (0.014  per- 
cent), manganese  (0.025  percent),  and  copper  (0.11  percent)  were 
found.  It  was  assumed  that  all  samples  maintained  these  purity 
limits  even  though  no  other  measurements  were  made. 
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I used  the  water  displacement  method  to  measure  an  ini- 
tial iron  density  of  7.868  + 0.012  gm/cm^.  Hardness  values 
ranged  from  69  to  71  on  the  Rockwell  F scale  for  two  differ- 
ent iron  samples  with  faces  surface-ground.  No  difference  in 
hardness  beyond  experimental  scatter  was  detected  along  diameter 
paths  on  the  faces  of  the  samples. 

An  unshocked  iron  sample  from  the  center  of  the  bar  was 
polished  and  etched  with  a 2 percent  Nital  solution  (2  cc  of 
concentrated  HNO^  and  100  cc  of  amyl  alcohol) . Grains  were 
irregular  in  shape  and  had  an  average  diameter  of  approximately 
0.1  mm.  Signs  of  foreign  material  were  evident  in  the  photo- 
micrograph. 

3.2.2.  Sample  Preparations 

Samples  were  cut  with  faces  perpendicular  to  the  bar 
axis.  Each  disc  was  turned  to  the  desired  diameter  and  its  faces 
were  surface-ground  to  the  desired  thickness,  which  ranged  from 
1 to  25.4  mm.  Faces  were  then  lapped  to  a flatness  of  less 
than  2 fringes  measured  with  sodium  light  and  an  optical  flat. 

An  exception  was  the  25. 4-mm-thick  plate  used  to  obtain  the 
equilibrium  transformation  stress.  This  iron  plate  would  not  fit 
the  available  lapping  wheel  and  was  therefore  hand  lapped  on 
600  grit  paper  placed  on  a flat  layout  table.  Its  final  flat- 
ness measured  less  than  20  fringes  using  sodium  light  and  an 
optical  flat.  The  two  faces  were  parallel  to  better  than  12  m 
for  all  samples  as  measured  by  a dial  indicator. 
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3.3.  Types  of  Experiments 

The  three  types  of  experiments  performed  in  this  work  are 
summarized  in  Table  3.1.  In  elastic  precursor  experiments, 
the  iron  samples  were  impacted  by  an  aluminum  projectile.  Stress 
histories  were  obtained  by  recording  with  oscilloscopes  the 
current  output  of  a quartz  gage  affixed  to  the  rear  face  of  the 
sample.  In  plate  slap  experiments  the  iron  samples  were  impacted 
by  an  aluminum  flier  plate.  In  the  explosively  driven  experi- 
ment, the  25 . 4-mm-thick  iron  sample  was  impacted  by  a detona- 
tion wave.  In  plate  slap  and  explosively  driven  experiments, 
shock  wave  transmission  times  and  free  surface  velocities  were 
determined  by  recording  with  a streak  camera  changes  in  the 
amount  of  light  being  reflected  from  the  surfaces  of  mirrors. 

All  experiments  were  designed  to  be  one-dimensional. 
Lateral  relief  from  the  sides  of  samples  during  the  recording 
time  was  avoided  by  maintaining  a sample  ratio  of  diameter  to 
thickness  greater  than  three.  The  only  two-dimensional  hydrody- 
namic effects  present  were  those  induced  by  tilted  or  unavoidably 
curved  shock  fronts. 

3.4.  Elastic  Precursor 
Experiments 

Two  elastic  precursor  experiments  were  performed  using 
the  gas  gun  facility  at  Washington  State  University.  Experi- 
mental details  for  these  were  essentially  the  same  as  reported 
41 

elsewhere. 
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Initiated  by  Pentolite  plane  wave  boosters. 
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3.5.  Plate  Slap  Experiments 

Shock  wave  transmission  experiments  were  performed  on 
thin  samples  of  Armco  iron  in  which  final  driving  stress  was 
near  200  kbar.  Specific  details  of  plate  slap  experiments  are 
given  in  the  following  subsections.  Included  are  discussions  of 
flier  system,  target  construction,  and  laboratory  methods. 


3.5.1.  Flier  Plate  System  for 
Plate  Slap  Experiments 

The  shock  wave  system  is  illustrated  in  Fig.  3.1.  The 
system  consists  of  a plane  wave  explosive  booster  (a) , a cast 
2,4,6-Trinitrotoluene  (TNT)  pad  (b) , a brass  plate  (c) , an  alumi- 
num flier  plate  (d) , a vacuum  chamber  (f,  g,  and  k) , and  a 
target  (g)  which  holds  the  iron  sample.  Reproducibility  of  the 
system  was  about  18  percent,  as  indicated  by  the  range  of  free 
surface  velocities  (1.01-1.19  mm/psec)  obtained  for  the  plas- 
tic II  shock  in  iron. 

The  difference  in  impedance  between  the  2 . 54-cm-thick 
brass  plate  and  the  0 . 32-cm-thick  aluminum  plate  causes  the 
brass-aluminum  interface  to  go  to  zero  pressure  when  backward- 
facing relief  waves  from  the  free  surface  of  the  aluminum  plate 
reach  the  interface.  The  aluminum  flier  plate  separates  from 
the  brass  plate  and  flies  free.  It  travels  12.7  mm  through 
vacuum  at  an  average  velocity  of  1.6  mm/psec  before  impacting 
the  iron  samples.  The  time  required  for  a flier  plate  to  travel 
the  12.7  mm  was  approximately  8 psec,  which  is  sufficient 
time  for  eight  wave  reverberations  to  occur  in  the  plate  before 
impacting  the  iron  sample. 
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Spalling  of  the  flier  plate  was  believed  not  to  occur 
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though  this  was  never  proved.  A graphical  solution  , ignoring 
Taylor  relief  waves  from  the  explosive  gas  products,  indicates  a 
maximum  tension  of  about  35  kbar  is  first  produced  near  the 
rear  boundary  of  the  flier  plate.  Therefore,  if  the  plate 
spalled  near  the  rear  surface,  it  would  not  affect  the  present 
results  because  the  thicker  front  of  the  plate  drives  the  shock 
for  a time  ample  for  all  measurements. 

A graphical  solution  treating  the  elastic-plastic  behav- 
ior of  aluminum  was  carried  out  for  times  t < d/c,  where  d is 
flier  plate  thickness,  c is  longitudinal  wave  velocity  in 
aluminum,  and  t=0  is  the  instant  of  separation  of  aluminum 
flier  from  the  brass  driver.  The  magnitude  of  the  tension  waves 
was  reduced  by  at  least  6 kbar  in  this  time.  The  attenuation 
occurs  because  elastic  waves  overtake  the  slower  plastic  relief 

waves.  The  physics  of  this  attenuating  process  has  been  treated 
4 3 

elsewhere.  At  the  time  of  impact  with  the  iron  sample  the 
aluminum  flier  plate  is  assumed  to  be  stress  free. 

_3 

All  flier  plates  were  tilted  less  than  4 x 10  rad 
relative  to  the  target  face,  in  a special  experiment  wherein  a 
transparent  glass  plate  was  used  in  place  of  the  iron  sample, 
flier  plate  impact  on  the  glass  plate  was  simultaneous  to 
0.012  (jsec  over  a central  area  with  diameter  of  30  mm.  Flier 
plate  edges  arrived  at  the  glass  plate  slightly  ahead  of  the 
center  portion,  but  not  enough  for  their  impact  to  perturb  the 
iron  samples  and  interfere  with  free  surface  measurements. 
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3.5.2.  Target  Preparation  for 
Plate  Slap  Experiments 

As  shown  in  Fig.  3.1,  the  target  consisted  of  a 15.24-cm- 
diameter  aluminum  ring  (g) , a 5.08-cm-diameter  aluminum  ring 
(j),  the  iron  sample  (e) , mirrors  (i),  and  epoxy  (h) . The 
15 . 24-cm-diameter  aluminum  ring  (target  ring  g)  served  to 
house  the  rest  of  the  components  and  as  part  of  the  vacuum  cham- 
ber. A 0 . 64-cm-thick  piece  of  plate  glass  which  served  as  the 
top  of  the  vacuum  chamber  allowed  mirrors  to  be  viewed  by  the 
streak  camera.  The  target  ring  was  machined  so  that  there  was  a 
small  angle  between  the  top  glass  and  the  iron  surface  to  insure 
that  light  reflections  from  glass  surfaces  would  not  reach  the 
camera.  The  target  ring  also  had  inserted  through  its  wall  a 
copper  tube  with  1.27-cm  inner  diameter  to  allow  pumping  of  a 
vacuum.  Pressures  of  less  than  0.03  bar  were  obtained. 

To  insure  reproducibility  of  the  shock  system,  all  dimen- 
sions of  the  component  parts  were  maintained  from  shot  to  shot 
within  0.005  cm.  Brass  plate,  aluminum  flier  plate,  and  spacing 
ring  surfaces  were  parallel  to  within  0.0038  cm.  The  brass 
plate  and  aluminum  flier  plate  faces  were  flatter  than  2 light 
fringes  as  measured  with  sodium  light  using  an  optical  flat. 
Mirrors  were  cut  from  a large,  optical  quality,  front-silvered 
mirror,  flat  to  1 light  fringe  as  measured  using  sodium  light 
and  an  optical  flat. 

Tilted  mirrors  were  placed  on  the  iron  samples  with  one 
edge  of  the  silvered  side  in  contact  with  the  iron  surface.  A 
contact  where  no  transmitted  light  was  visible  was  considered 
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satisfactory.  The  contact  edge  of  the  mirror  was  placed  far 
enough  away  from  the  edge  of  the  iron  sample  to  avoid  interfer- 


ence from  edge  effects.  The  opposite  end  of  the  mirror  was 
propped  up  by  a jeweler's  drill  to  obtain  the  desired  angle 
(typically  2 deg)  with  the  iron  surface.  The  jeweler's  drill 
was  removed  after  the  mirror  was  fixed  securely  in  place  with 
epoxy.  Mirrors  with  small  angles  of  inclination  were  used  to 
insure  that  the  gap  between  iron  and  mirror  surfaces  closed  be- 
fore wave  interactions  in  the  mirrors  could  affect  the  measure- 
ments. 

The  angle  between  tilted  mirror  and  iron  su.face  was 
measured  with  a shadowgraph  instrument  which  is  accurate  to 
1 min  of  arc.  The  iron  impact  surface  opposite  the  surface 
with  mirrors  was  used  as  a reference  plane  for  the  angle  measure- 
ment. Since  iron  surfaces  were  parallel  to  better  than  12 
any  errors  in  the  measured  angle  due  to  this  lack  of  perfect 
parallelism  were  smaller  than  could  be  measured.  Some  tilted 
mirror  angles  were  also  measured,  using  a spectrometer  table, 
by  observing  with  a telescope  the  super imposi t ion  of  a cross-hair 
and  its  image  reflected  alternately  from  the  tilted  mirror  and 
the  flat  mirror  against  the  iron  surface.  Angles  could  be  mea- 
sured in  this  way  to  within  0.5  min  of  arc.  The  remeasured 
angles  agreed  within  1.5  min  of  arc  with  shadowgraph  measure- 
ments. 

The  target  assembly  of  Fig.  3.1  was  similar  to  that  for 

4 1 

precursor  experiments  discussed  elsewhere.  A target  blank  was 
built  by  wringing  the  target  ring  (g)  down  against  a flat  plate. 
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The  flat  plate  was  hand  lapped  on  600  grit  paper  to  a flatness 
better  than  2 sodium  light  fringes  as  measured  with  an  optical 
flat.  The  small  ring  (j)  was  also  wrung  down  against  the  flat 
plate  near  the  center  of  the  larger  ring.  The  small  ring  even- 
tually housed  the  iron  sample  (e) . Vacuum  grease  between  sur- 
faces and  on  the  entire  face  of  the  flat  plate  insured  easy 
removal  of  the  ring  from  the  plate  after  the  epoxy  had  been 
poured  and  cured.  Epoxy  (Shell  Epon  815  mixed  4/1  by  weight  with 
a T-1  hardener)  was  poured  inside  the  large  ting  to  a thickness 
of  1.27  cm.  The  smaller  ring  was  left  blank  at  this  stage  of 
target  preparation.  The  epoxy  was  poured  in  increments  of  50  gm 
and  each  increment  was  allowed  to  cure  before  the  next  was 
poured.  The  target  blank  was  then  removed  from  the  flat  plate 
and  vacuum  holes  were  drilled  in  the  epoxy  to  allow  air  to  flow 
to  the  vacuum  line  from  the  chamber  where  the  flier  plate  was 
located. 

The  target  blank  was  hand  lapped  on  600  grit  paper  on 
a layout  table  to  obtain  a flat  surface;  it  was  then  wrung  down 
on  the  flat  plate  with  vacuum  grease.  Two  0.64-cra-wide  and 
1.27-cm-long  mirrors  were  wrung  down  against  the  flat  steel 
plate  next?  to  the  iron  sample  to  measure  the  flier  arrival  and 
flier  tilt.  The  iron  sample  and  flat  mirrors  in  the  impact  plane 
were  then  epoxied  in  place.  The  first  three  or  four  experiments 
used  the  Shell  epoxy  to  fix  the  sample  and  mirrors  in  place.  It 
was  discovered  that  for  thin  samples  (1.5  mm  or  less)  a 
stronger,  more  viscous  epoxy  (Epoxi-Patch  Kit  1C)  was  better 
suited  to  this  purpose.  The  completed  target  was  then  removed 
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from  the  flat  plate.  Further  lapping  could  not  be  done  at  this 
point  because  the  silver  would  be  removed  from  the  mirror  sur- 
faces located  at  the  impact  plane. 

This  lack  of  a final  lapping  contributed  slightly  to  the 
tilt,  since  sample  and  mirrors  could  not  be  put  down  perfectly 
flat.  Indeed,  tilts  observed  from  this  flier  system  were  1 
order  of  magnitude  greater  than  that  obtained  in  the  precursor 
experiments.  Not  all  of  this  can  be  ascribed  to  lack  of  the 
final  lapping;  there  are  many  factors  other  than  this  contribut- 
ing to  tilt.  I concluded  that  imperfections  in  design  and  con- 
struction of  the  explosive  plane  wave  booster  were  the  main 
source  of  the  observed  tilt.  These  were  beyond  the  control  of 
the  author,  and  it  was  therefore  concluded  that  tolerances  and 
methods  of  target  preparation  used  were  adequate.  A better  plane 
wave  booster  would  be  required  to  justify  closer  tolerances  in 
target  preparation. 

3.5.3.  Laboratory  Methods 

The  target  and  explosive  charge  were  mounted  inside  the 
test  chamber  on  a leveled  table.  The  table  was  placed  so  that 
the  target  assembly  could  be  viewed  through  the  slit  of  a rotat- 
ing mirror  streak  camera  aligned  along  the  centers  of  the  tilted 
mirrors.  A diagram  of  the  setup  is  shown  in  Fig.  3.2.  Light 
from  an  exploding  wire  light  source  was  reflected  from  mirror 
surfaces  located  in  the  target  and  intercepted  by  the  camera. 
Incident  light  on  the  mirrors  was  at  a slight  angle  from  the 
normal  to  the  target  to  allow  convenient  placement  of  the  light 
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source.  Alignment  of  the  experiment  with  the  camera  was  accom- 
plished by  projecting  light  through  the  camera  lens  and  slit  out 
into  the  test  chamber,  where  it  was  incident  on  the  reflecting 
surfaces  of  the  mirrors.  These  surfaces  then  reflected  the 
light,  at  a small  angle  from  the  normal  to  the  target  surface, 
onto  the  glass  capillary  tube  surface  of  the  light  source.  This 
insured  that  light  source  and  mirrors  were  properly  aligned  with 
the  camera. 

The  Cordin  132  streak  camera  was  outside  the  explosive 
chamber  and  protected  from  blast  by  a 2. 54-cm-thick  glass  view- 
ing window.  The  camera  was  focused  on  a precision  grid  placed 
on  the  surface  of  the  target.  A still  picture  was  taken  and 
developed  before  every  experiment  to  verify  focus.  The  camera 
was  also  run  at  the  same  writing  speed  used  in  the  experiment 
and  an  aligned  light  wire  source  was  exploded  to  check  the  amount 
of  light  available  for  the  experiment. 

The  exploding  wire  light  source  was  a 0 . 005-cm-diameter 
tungsten  wire  with  resistance  of  3.4  ohm,  inside  a 10.15-cm- 
long  glass  capillary  tube  having  an  outer  diameter  of  0.74  cm 
and  an  inner  diameter  of  0.15  cm.  Experience  has  shown  this  to 
be  a good  light  source.  Air  inside  the  capillary  was  shocked  to 
luminescence  by  discharging  a 3 uF  capacitor,  which  has  been 
charged  to  8000  V,  through  the  tungsten  wire.  The  capacitor 
power  supply  was  connected  to  the  wire  by  a 3.7-m-long, 

10,000-V  test  prod  wire  (18  AWG) . 

Kodak  TRI-X  70-mm-wide  film  was  used  in  the  streak  cam- 
era. It  was  developed  in  Kodak  D-11  developer  and  fixed  in  Kodak 
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Rapid  fixer.  The  film's  dynamic  spatial  resolution  for  this 
experimental  setup  of  >300  lines/cm  was  measured  in  a special 
experiment  by  successfully  recording  the  image  of  a Ronchi  ruling 
as  a function  of  time  for  a camera  writing  rate  of  15  mm/usec. 
Time  resolution  of  the  camera  is  a function  of  camera  speed, 
rotor  mirror  distortion,  spatial  resolution,  slit  width,  film 
characteristics,  and  ability  to  read  the  film  record.  A good 
estimate  of  time  resolution  is  obtained  by  dividing  slit  width 
by  camera  speed.  This  estimate  gives  a time  resolution  of 
0.0033  lisec  for  a camera  writing  speed  of  15  mm/psec  and  a 
slit  width  of  0.05  mm.  Possible  errors  in  data  due  to  time 
resolution,  spatial  resolution,  and  other  such  reasons  are  dis- 
cussed in  Appendix  C. 

3.6.  Explosively  Driven 
Experiments 

An  explosive  shock  system  in  contact  with  a 25.4-mm- 
thick  iron  sample  was  used  to  obtain  free  surface  data  for  driv- 
ing stress  near  250  kbar.  The  explosive  shock  system  consisted 
of  a plane  wave  booster  in  contact  with  a 2. 54-cm-thick  pressed 
ammonium  perchlorate  (AP)  explosive  pad.  The  AP  particle  diam- 
eters were  less  than  10  p before  pressing  to  a density  of 
1.3  gm/cm^.  Paces  of  the  charge  were  machined  to  parallelism 
better  than  0.0036  cm  over  a 12.7-cm  diameter.  Paces  of  the 
charges  were  gently  hand  lapped  on  a layout  table  with  600  grit 

paper  to  insure  a flat  surface. 

An  AP  pad  was  chosen  principally  because  of  the  250  kbar 

stress  it  induces  in  iron  and  the  observed  improvements  in 
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planarity  of  the  shock  front,  compared  to  the  shock  from  the  !] 

plane  wave  booster  along.  For  example,  arrival  of  a shock  wave  I 

at  the  free  surface  of  a 1. 27-cm-thick  aluminum  plate  driven 
by  this  explosive  system  was  simultaneous  to  better  than 

0.04  ysec  over  an  area  with  a 5.93-cra  diameter.  Typical  | 

shock  wave  tilts  resulting  from  this  explosive  system  were 
5 X 10  ^ rad. 

Laboratory  methods  and  target  preparations  for  these 
experiments  were  similar  to  those  for  plate  slap  experiments. 
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CHAPTER  4 

EXPERIMENTAL  RESULTS 

Experimental  results  for  two  precursor  experiments,  one 
transmission  experiment  on  a thick  sample,  and  eight  transmis- 
sion experiments  on  thin  samples  are  reported  in  this  chapter. 
Stresses  behind  shock  fronts  were  inferred  from  measured  free 
surface  velocity  and  average  shock  velocity.  Elastic  data 
required  for  calculation  of  stress  behind  the  plastic  I shocks 
were  taken  from  Taylor  and  Rice's  paper  on  Armco  iron.^^ 

Errors  have  been  reported  in  the  shorthand  form  "A  + e" 
where  e is  the  average  deviation  of  A.  Appendix  C treats  the 
error  analysis  in  some  detail. 

4.1.  Elastic  Precursor  Data 

Two  precursor  experiments  were  performed  for  different 
driving  stresses.  Sample  thicknesses  were  3.1  and  6.3  mm; 
results  are  summarized  in  Table  4.1  and  in  Fig.  4.1.  The  final 
stresses  in  iron  were  determined  by  the  intersection  of  the 
equilibrium  P-u  curve  of  iron  and  the  reflected  P-u  curve  of 
aluminum,  the  aluminum  curve  being  centered  at  P=0  and  u = 
projectile  velocity. 

The  solid  line  in  Fig.  4.1  was  inferred  from  Taylor  and 
Rice's^^  free  surface  velocities.  Precursor  stresses  reported 
here  are  slightly  lower  than  the  solid  line  in  Fig.  4.1. 
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FIGURE  4.1  ELASTIC  PRECURSOR  STRESS  AS  A FUNCTION  OF  SAMPLE  THICKNESS. 
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Deviations  of  these  two  data  points  from  the  line  are  within 
experimental  error;  discussion  will  be  deferred  to  Section  6.2. 

4.2.  Transmission  Data  for 
a 25. 4-mm-Thick  Sample 

A single  experiment  was  performed  to  obtain  stress 
behind  the  plastic  I shock  on  a thick  sample.  Results  are  used 
as  a reference  for  thin  sample  results.  A stress  of  131.4  + 

3.3  kbar  was  inferred  from  measured  free  surface  velocity  of 

0.653  + 0.0117  mm/ysec  and  an  average  measured  shock  velocity 

of  5.082  + 0.090  mm/usec.  A Hugoniot  elastic  limit  (HEL)  of 

36 

8.3  kbar  was  used  for  the  plastic  I stress  calculation.  This 

2 

result  agrees  well  with  Bancroft,  Peterson,  and  Minshall's 
value  of  131.0  kbar  for  a 24 . 7-mm-thick  Armco  iron  sample. 

4.3.  Transmission  Data  for 
Thin  Samples 

Eight  transmission  experiments  were  performed  at  approxi- 
mately 201-kbar  driving  stress.  Sample  thicknesses  ranged  from 
0.94  to  6.31  mm.  Free  surface  motion  and  average  shock  veloci- 
ties for  the  iron  samples  were  measured. 

4.3.1.  Plastic  I Data 

Table  4.2  summarizes  results  for  the  plastic  1 shock. 

The  table  includes  sample  thicknesses,  free  surface  velocities, 
average  shock  velocities,  and  shock  wave  transit  times.  Depend- 
ence of  free  surface  velocity,  stress,  and  stress  jump,  P^  - 
on  sample  thickness  is  shown  in  Figs.  4.2  to  4.4.  These  plots 
show  small  dependence  of  free  surface  velocity  and  stress  on 
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FIGURE  4.2  FREE  SURFACE  VELOCITY  OF  IRON  AFTER  PLASTIC  I SHOCK  ARRIVAL. 
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sample  thickness.  The  plot  of  stress  jump  versus  thickness 
implies  that  the  increase  in  stress  behind  the  plastic  I front 
for  small  sample  thicknesses  may  be  accounted  for  by  the  known 
increase  in  stress  of  the  precursor  as  sample  thickness  is 
decreased. 

Figure  4.2  shows  free  surface  velocity  versus  sample 
thickness;  data  from  two  additional  sources  are  shown  for  com- 
parison. These  other  data  points  result  from  experiments  using 
different  experimental  techniques  and  different  final  driving 
stresses.  The  solid  line  is  a least-squares  fit  of  present  data 
to  a straight  line.  It  corresponds  to  an  increase  in  free  sur- 
face velocity  of  about  0.1  percent  for  every  1-mm  decrease 
in  sample  thickness.  Its  slope  is  not  significantly  different 
from  zero. 

Stresses  behind  the  plastic  I shock  in  iron  versus  sample 

thickness  are  shown  in  Fig.  4.3.  In  addition,  two  solid  curves 

are  present  which  are  calculated  from  the  phenomenological  model 

20 

of  Duvall  and  Hone  for  two  different  relaxation  times.  Com- 
parison between  calculated  curves  and  experimental  data  shows 
that  relaxation  time,  x,  is  obviously  less  than  0.1  ysec. 

The  significance  of  this  will  be  discussed  in  Section  6.1. 

The  data  in  Fig.  4.3  show  a slight  increase  in  stress 
for  decreasing  sample  thickness.  Stress  is  about  140  kbar  for 
a sample  thickness  of  1 mm  and  about  136  kbar  at  6.35  mm. 

The  140-kbar  value  is  9 kbar  greater  than  the  stress  meas- 
ured for  the  25.4-mm-thick  sample  reported  in  Section  4.2. 
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This  increase  in  stress  for  decreasing  sample  thickness  is  not 
very  significant. 

Stress  jumps  across  the  plastic  I shock  versus  sample 
thickness  are  shown  in  Fig.  4.4.  Precursor  amplitudes  were 
taken  from  Fig.  4.1.  The  solid  line  represents  the  weighted 
average  of  the  stress  jump  for  the  present  results.  The  stress 
jump  appears  constant.  Deviations  from  this  constant  value  of 
123.6  + 1.6  kbar  for  all  the  data  are  within  experimental  error. 
Constancy  of  the  stress  jump  implies  that  the  increase  in  plas- 
tic I stress  for  decreasing  sample  thickness  shown  in  Fig.  4.3 
may  be  due  to  the  increase  in  precursor  stress. 

4.3.2.  Plastic  II  Data 

Table  4.3  summarizes  results  for  the  plastic  II  shock  in 
iron.  These  results  are  from  the  same  experiments  as  those  that 
produced  the  plastic  I data  of  Table  4.2.  The  table  includes 
sample  thicknesses,  free  surface  velocities,  shock  wave  trans- 
mission times,  and  rise  times  of  the  plastic  II  shock  front. 
"Beginning"  of  the  plastic  II  wave  front  refers  to  the  time  at 
which  free  surface  velocity  begins  to  exceed  the  plastic  I value. 
"Top  of  the  wave"  refers  to  the  time  when  free  surface  velocity 
again  becomes  constant  or  nearly  so.  Plastic  II  results  are 
used  here  primarily  to  determine  a best  value  for  Lagrangian 
velocity  of  the  wave  front  and  to  obtain  the  best  possible 
pressure-volume  (P-V)  point  from  the  data  set. 

Plastic  I and  plastic  II  shock  wave  arrival  times  at  the 
free  surface  for  various  sample  thicknesses  are  shown  in  Fig.  4.5, 
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The  arrival  times  were  fit  by  least  squares  to  straight  lines, 
with  arrival  time  taken  as  the  independent  variable  and  sample 
thickness  as  the  dependent  variable.  Each  arrival  time  was 
assigned  a weighting  factor  proportional  to  the  inverse  square 
of  its  error.  Differences  in  final  stress  reached  in  the  iron 
samples  were  ignored  when  determining  Lagrangian  velocity  of  the 
plastic  II  shock  wave. 

Lagrangian  velocities  obtained  by  this  procedure  are; 

I 

plastic  I wave,  U2  = 5.074  + 0.045  mm/psec 

t 

beginning  of  plastic  II  wave,  U^g  = 3.883  + 0.138  mm/ysec; 

I 

top  of  plastic  II  wave,  “ 4.121  + 0.235  mm/ysec 

The  beginning  plastic  II  wave  path  intersects  the  h=0  plane  at 
t = 0.23  + 0.017ysec.  The  top  of  the  wave  front  appears  to 
propagate  slightly  faster  than  the  bottom,  but  errors  in  the 
velocities  do  not  allow  any  significance  to  be  assigned  to  this 
observation.  In  fact,  the  error  bars  are  of  such  large  magni- 
tude that  no  conclusion  about  steadiness  of  the  wave  front  can 
be  drawn  from  the  data  of  Fig.  4.5. 

The  h-t  plots  of  shock  transmission  times  in  Fig.  4.5 
do  not  represent  the  true  paths  of  the  initial  forward-facing 
plastic  shocks  because  of  interactions  with  backward-facing 
waves  from  the  free  surfaces,  described  in  Section  2.2.  Correc- 
tions to  the  plastic  I path  are  negligible  because  of  the  small 
amplitude  of  the  precursor  reflection.  Values  of  (h,t)  points 
for  the  beginning  of  the  plastic  II  wave,  corrected  according 
to  Eqs.  (2.9)  and  (2.10),  are  shown  in  Fig.  4.6.  The  rectangles 
represent  estimated  errors  for  each  point.  The  fitted  line 
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FIGURE  4.6  CORRECTED  LAGRAGIAN  h-t  DIAGRAM  FOR  BEGINNING  OF  THE  PLASTIC  II  FRONT. 
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corresponds  to  a velocity  of  3.587  + 0.227  nun/ysec,  which  is 
slightly  less  than  that  obtained  for  the  uncorrected  data  of 
Fig.  4.5. 


4.3.3.  Stress-Volume  State 
Behind  the  Plastic  II  Shock 

The  stress-volume  state  behind  the  plastic  II  shock  was 
calculated  using  the  equilibrium  jump  equations,  Eqs.  (2.4)  and 
(2.5);  weighted  averages  of  plastic  I and  plastic  II  particle 
velocities;  and  Lagrangian  shock  velocities  from  the  preceding 
section ; 

Uj^  = 0.032  + 0.005  mm/ysec, 

U2  = 0.333  + 0.004  mm/ysec, 

Uj  = 0.565  + 0.018  mm/ysec, 

I 

* 6.18  + 0.21  mm/ysec, 

I 

U2  = 5.074  + 0.045  mm/ysec, 

= 3.587  + 0.227  mm/ysec, 

P3  = 201  + 8.4  kbar, 

V3/V0  = 0.871  + 0.008, 

Vq  = 7.098  + 0.011  cc/mol. 

This  stress-volume  point  is  given  in  Fig.  4.7  along  with  data 

2 15 

from  Bancroft,  e^  data  from  Barker  and  Hollenbach,  and 

27 

Andrews'  equilibrium  curve.  It  is  reasonably  consistent  with 
the  Bancroft  data  and  the  Barker  data. 


I-, 

'5 
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4.4.  Summary 

The  most  significant  experimental  result  of  this  study 
is  shown  in  Fig.  4.3  and  4.4:  There  is  little  or  no  variation 
of  plastic  I wave  amplitude  for  propagation  distances  between 

0.9  and  6.35  mm.  This  implies,  according  to  Fig.  4.3,  a 
relaxation  time  of  about  0.05  ysec  foe  onset  of  the  a -►  e 
transition.  When  the  single  measurement  of  25.4-mm  distance 
is  included  with  close-in  data,  one  infers  a slow  variation  of 
transition  stress  with  distance,  superimposed  on  the  rapid  decay 
below  1 mm.  Inference  of  the  initial  decay  depends  on  the 
assumption  that  initial  compression  at  the  impact  surface  was 
entirely  in  the  a phase. 

Further  results,  which  are  essentially  corrobative,  are 

that : 

1.  Elastic  precursor  amplitude  increases  as  sample 
thickness  is  decreased. 

2.  Transition  stress  measured  in  a 25 . 4-mm-thick 
sample  is  131.4  + 3.3  kbar . 

3.  Relative  volume  behind  the  plastic  II  wave  is 

V3/V0  = 0.871  + 0.008  for  a stress  of  201  + 8.4  kbar. 

4.  Rise  time  in  the  plastic  II  shock  front  is 
0.18  + 0.02/isec. 

These  results  are  in  essential  agreement  with  those  of  earlier 
experiments  reported  in  references  36,  2,  2,  and  23,  respectively 
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CHAPTER  5 

THEORY  OF  THE  ALPHA  TO  EPSILON  TRANSFORMATION 

In  this  chapter  shock  wave  measurements  on  iron  and 
metallurgical  data  on  martensitic  transformations  are  brought 
together  in  an  attempt  to  produce  a plausible  description  of  the 
iron  transformation. 

5.1.  Review  of  Experimental 
Informat  ion 

The  alpha  to  epsilon  transformation  in  iron  has  been 
studied  by  several  different  investigators  using  different  ex- 
perimental techniques  and  varying  different  parameters.  Salient 
results  of  these  studies  are  summarized  here. 

5.1.1.  Shock  Experiments 

Stress  in  the  plastic  I shock  is  a measure  of  the  trans- 
TL 

formation  stress,  P , for  the  compression  or  loading  process. 

TL 

Some  measurements  of  P for  19-  to  25-mm-thick  iron  samples 

are  given  in  Table  5.1.  Iron  samples  used  have  had  different 

metallurgical  histories  and  probably  contained  different  trace 

amounts  of  impurities.  Agreement  among  values  of  P in 

Table  5.1  implies  that  trace  amounts  of  impurities  and  different 

TL 

metallurgical  histories  do  not  affect  P significantly. 

22 

Minshall  has  shown  that  work  hardening  and  heat  treating  Armco 

TL 

iron  slightly  increases  P , but  the  transition  volume  remained 
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constant.  Although  trace  amounts  of  impurities  have  little 

TL  TL 

effect  on  P , 1.5  percent  of  carbon  increases  P by 

44 

5 kbar.^^ 

Barker  and  Hollenbach^^  measured  shock  release  states  in 

iron  and  showed  that  the  reverse  transformation  of  epsilon  to 

TU 

alpha  iron  is  initiated  at  a stress,  P , of  98  kbar.  The 

TL 

stress  for  equilibrium  transformation  lies  between  P and 

Til  T n 

P . -Their  mean  value  is  P = 116  kbar,  for  which  T = 332  K 
on  the  Hugoniot. 

TL 

The  6 percent  decay  in  P for  wave  propagation  from 

6 to  25  mm  described  in  Chapter  4 is  similar  to  other  reported 

2 15 

values.  ' Measurements  on  specimens  down  to  a thickness  of 
1 mm  reported  in  Chapter  4 show  that  initial  decay  of  the  plas- 
tic I shock  stress  is  very  rapid.  This  observation  is  supported 
by  Barker  and  Hollenbach  in  measurements  on  specimens  of  thick- 
nesses of  3 mm  and  greater. 

Shock  compression  data  on  iron  show  that,  for  a given 

volume,  measured  stress  in  the  mixed  phase  region  exceeds  that 

27  28 

expected  for  the  equilibrium  Hugoniot.  ' (See  Fig.  4.7.) 
Transformation  begins  at  130  kbar  and  appears  to  be  complete 

near  200  kbar.  The  slope,  -5  Mbar  gm/cm  (1  Mbar  = 

3 TL 

10  kbar),  of  the  Hugoniot  just  above  P disagrees  signifi- 
cantly with  the  slope  for  the  equilibrium  curve,  which  is  about 
-0.1  Mbar  gm/cm^.  The  equilibrium  value  is  calculated  for 
conditions  of  uniform  hydrostatic  pressure  throughout  the 
material.  For  various  reasons,  including  anisotropy  of  stress, 
measured  values  of  stress,  P,  differ  from  hydrostatic  values. 


58 


NSWC/WOL  TR  77-137 


P,  by  unknown  amounts,  which  may  account  for  some  of  the  differ- 
ences in  slopes.  It  is  not  certain  that  measured  Hugoniot  states 
in  the  mixed  phase  region  are  in  equilibrium.  They  persist  for 
the  experimental  time  span  and  will  be  denoted  "quasistatic" 
states. 

Values  of  P for  various  temperatures  are  given  in 

Fig.  5.1.  The  Clausius-Clapeyron  equilibrium  relation  and  the 

equation  of  state  given  in  Appendix  A yield  a value  of 

-60.8  X 10~®  Mbar/°K  for  ^ . This  calculated  slope  is  within 

45 

experimental  error  of  the  observations  of  Johnson,  et  al. 

5.1.2.  Static  Experiments 

A number  of  static  isothermal  compression  experiments  on 

iron  have  shown  that  the  transformation  is  initiated  near 

130  kbar  and  goes  to  completion  for  stresses  greater  than 

145  kbar Giles,  et  al.^^  found  values  of  133  kbar 

for  p"^^  and  81  kbar  for  p"^^  along  a 300°K  isotherm. 

They  also  found  that  the  alpha  phase  persisted  for  stresses 

above  163  kbar  while  Mao,  et  al. report  no  signs  of  alpha 

32 

iron  above  145  kbar.  Giles,  ^ al.  found  that  the  epsilon 
phase  persisted  down  to  45  kbar  in  the  unloading  process.  The 
mean  value  for  transformation  stress,  reported  by  Giles, 

et  al.  was  107  + 8 kbar,  which  is  9 kbar  less  than  the  mean 
value,  reported  by  Barker  and  Hollenbach . Some  of 

this  difference  may  be  due  to  different  amounts  of  shear  stress 
in  static  and  dynamic  experiments. 
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FIGURE  5.1  TRANSFORMATION  STRESS  FOR  IRON  UNDER  COMPRESSION 
AT  VARIOUS  TEMPERATURES. 


60 


j 


1 


NSWC/WOL  TR  77-137 


I 


5.1.3.  Discussion 

Differences  between  the  measured  Hugoniot  and  calculated 
equilibrium  Hugoniot  may  result  from  several  causes  including: 

(1)  partial  transformation,  (2)  transformation  to  a metastable 
state,  (3)  surface  energy  requirements  for  the  new  phase  to 
nucleate,  (4)  strain  energy  requirements  for  accommodation  of 
the  new  phase  in  the  parent  lattice,  (5)  effects  of  shear  on  the 
transformation,  and  (6)  pre-existing  new  phase  nuclei  of  various 
sizes. 

The  difference  between  pTeq  presumably 

due  to  shear.  However,  there  exists  no  clear  evidence  that 
shear  in  iron  is  important  to  the  phase  transformation,  except 
for  effects  of  work  hardening  and  heat  treating  prior  to  shock 
loading . 

Several  investigators  have  suggested  that  the  alpha  to 

32 

epsilon  transformation  is  martensitic.  Giles,  et  al.  make  the 

—TL 

strongest  case  based  on:  (1)  the  large  difference  between  P 
and  P , (2)  existence  of  two  phases  over  a large  pressure 
range,  and  (3)  simultaneous  transformation  of  independent  areas. 
These  properties  ate  all  characteristic  of  martensitic  trans- 
formations. 

5.2.  Characteristics  of  Marten- 
sitic Transformation 

Martensitic  transformations  are  defined  in  terms  of 
process,  not  in  terms  of  structure  of  properties.  They  require 
no  diffusion  of  atoms;  they  are  displacive  or  shearlike,  in  that 
the  atoms  move  cooperatively  to  produce  substantial  shape  changes 
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in  the  transforming  region,  even  though  each  atom  shifts  only 
slightly  relative  to  its  neighbors.  These  diffusionless  and  dis- 


placive  characteristics  distinguish  martensitic  transformations 

from  all  other  types.  ’ 

Martensitic  tranformations  produce  lenticular  plates 
with  semicoherent  interfaces  between  plates  and  parent  lattice. 

Plates  grow  with  a velocity  near  the  speed  of  sound  but  stop 


growing  when  (1)  a grain  boundary,  (2)  another  martensitic  plate, 

I 

or  (3)  other  lattice  disturbances  which  serve  as  barriers  are  i 

I, 

encountered.  The  distribution  of  plate  sizes  within  a grain  is  j 

not  well  understood.  Fisher^®  suggests  that  plates  form  randomly  | 

throughout  each  grain,  constantly  subdividing  grains,  so  later  i 

generations  of  plates  form  in  smaller  and  smaller  parent  volumes. 

47  i 

Magee  has  proposed  a somewhat  different  order  of  growth  in 

which  nucleation  is  not  random  throughout  each  grain,  so  at  any 

time  some  grains  will  contain  martensitic  plates  and  some  will  i 

j 

not.  These  plates  tend  to  occur  in  clusters,  with  various-sized  j 

plates  in  each  cluster.  It  has  been  observed  that  the  first  i 

detectable  amount  of  transformation  is  due  to  the  presence  of 

47  ' 

single  clusters  of  plates  in  several  grains,  in  contrast  with 

ii 

Fisher's  model.  Further  transformation  largely  involves  spread-  , 

i' 

ing  of  clusters  to  untransformed  regions.  To  a first  approxima-  i 

tion,  volume  of  newly-formed  plates  in  a cluster  is  constant.  ' 

I ■ 

This  picture  leads  naturally  to  the  conclusion  that  increments  [i 

in  volume  of  the  new  phase  are  proportional  to  increments  in  the  !< 

number  of  nucleation  centers  activated.  I 
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5.2.1  Crystallography 

Martensitic  growth  causes  bulk  shape  changes,  from  such 


causes  as  tilting  of  the  martensitic  plates,  which  are  visible 
at  flat  surfaces.  The  plates  are  surrounded  by  untrans- 
formed matrix  and  are  usually  lenticular,  although  in  some 


steels  "needles"  or  "laths"  have  been  reported.  These  plates 
may  in  turn  contain  a fine  structure  of  slip  bands  or  twins  as 
a direct  consequence  of  the  transformation  process.^® 

The  martensitic  plate  has  mirror  symmetry  about  a plane, 
and  this  plane  has  a particular  and  reproducible  orientation 
with  respect  to  the  parent  phase  for  any  lattice.  The  plane  in 
the  parent  lattice  which  lies  parallel  to  the  symmetry  plane  of 
the  martensitic  plate  is  called  the  habit  plane.  It  is  the 
plane  along  which  the  principal  shear  displacement  occurs  in  the 
martensitic  reaction. 

Consideration  of  experimental  data  on  martensitic  trans- 
formations has  encouraged  development  of  crystallographic  expla- 
nations for  any  particular  transformation  of  orientation  rela- 
tionships, habit  planes,  deformation  shapes,  and  nature  of  the 
fine  structure  within  the  martensitic  plates. 

Bowden,  et  al . concluded  that  markings  from  recovered 
shocked  alpha  iron  samples  resulted  from  shear  transformation  to 
the  epsilon  phase,  body-centered-cubic  (bcc) , to  hexagonal-close- 
packed  (hep)  and  its  reversal,  which  suggests  that  the  transfor- 
mation is  martensitic.  The  shear  mechanism,  reviewed  below,  is 

49 

essentially  the  same  as  proposed  by  Burgers  for  the  zirconium 
transformation. 
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Some  comments  on  notation  are  required.  Lattice  planes  | 

I I * 50  1 

are  defined  by  Miller  indices  with  integer  components  (a  ,b  ,c  ).  j 

Directions  in  the  lattice  are  given  as  vectors  with  integer  com-  j| 

I ■’ 

ponents  |a,b,cj.  A bar  above  a component,  such  as  (a  ,b  ,c  ), 
indicates  a negative  value.  In  cubic  crystals,  a direction 

ja,b,c2  is  perpendicular  to  a plane  (a,b,c)  having  the  same  | 

indices.  j 

A bcc  cell  has  six  different  symmetric  planes  which 
are  "most  densely  packed"  (not  close-packed):  (lIO),  (110), 

(101),  (Oil),  (Oil),  and  (lOl) . Figure  5.2(a)  shows  a most 
densely  packed  (110)  plane  with  rows  of  close-packed  atoms 
such  as  A,  B,  and  C along  each  cube  diagonal.  It  is  along 
these  close-packed  rows  that  slip,  twinning,  and  some  phase 
transformations  take  place  due  to  applied  shearing  forces.  These 
rows  lie  in  either  (Il2)  or  (ll2)  planes;  e.g.,  the  row 

ABC  in  Fig.  5.2(a)  is  in  the  (ll2)  plane.  | 

Changes  in  a particular  triclinic  cell  illustrate  how  ^ 

the  bcc  lattice  can  be  transformed  to  hep  by  shear  deforma- 
tion. In  Fig.  5.2(b)  three  sequential  (110)  planes  are  repre- 
sented in  an  isometric  drawing;  separations  between  planes  are  | 

exaggerated  for  clarity.  These  planes  are,  for  example,  those 
containing  atoms  D,  B,  and  E,  respectively,  as  shown  in 
Fig.  5.2(a).  The  triclinic  cell  of  interest  is  outlined  by 
dashed  lines  in  Fig.  5.2(b)  with  four  sides  perpendicular  to  the 
(110)  planes  and  lying  in  either  the  (Il2)  or  (112)  planes. 

It  is  along  these  sides  that  shearing  action  for  the  transforma- 
tion will  take  place.  The  shear  mechanism  for  transformation  is 
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easy  to  visualize  if  the  triclinic  cell  is  oriented  as  shown  in 

Fig.  5.3(a).  To  obtain  this  orientation  of  the  cell  requires 

three  rotations  of  the  lattice  with  respect  to  the  set  of  fixed 

coordinates  in  Fig.  5.2(a):  (1)  rotate  the  lattice  around  the 

z axis  through  an  angle  of  +45°,  which  gives  the  parallel 

(110)  planes  shown  in  Fig.  5.2(b);  (2)  rotate  the  lattice 

around  the  y axis  by  +90°,  which  results  in  the  (110) 

planes  parallel  to  the  xy  plane;  (3)  rotate  the  lattice  around 

o * 

the  z axis  through  an  angle  of  -35  15  , which  orients  the 
triclinic  cell  as  in  Fig.  5.3(a). 

The  shear  mechanism  for  bcc  to  hep  transformation  is 
illustrated  in  Fig.  5.3(a)  and  (b)  where  a shear  on  the  (112) 
plane  in  the  [^lllj  direction  changes  the  70°30  angle  between 
the  two  sides  in  the  basal  plane  to  a 60°  angle.  The  open 
circles  in  Fig.  5.3(b)  represent  atomic  positions  before  shear. 
The  resulting  cell  has  hexagonal  symmetry  but  is  not  close-packed 
because  of  the  central  atom.  The  central  atom  must  be  shifted, 
as  indicated  by  the  two  open  circles  in  the  center  plane  of  the 
cell  of  Fig.  5.3(b),  to  bring  it  into  line  with  the  interstitial 
position  between  the  trio  of  atoms  in  the  basal  planes  above  and 
below.  Further  slight  adjustments  of  cell  dimensions  are  neces- 
sary for  conformance  with  lattice  parameters  of  the  hep  phase. 

The  transition  is  initiated  when  a critical  value  of 
shear  stress  is  exceeded  for  a row  of  close-packed  atoms  such  as 
the  direction  on  the  (112)  plane. 
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(il2)  PLANE 


FIGURE  5.3  SHEAR  MECHANISM  FOR  BODY-CENTEREO-CUBIC  TO 
HEXAGONAL  CLOSE-PACKED  TRANSFORMATION. 
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L. 


5.2.2.  Uniaxial  Strain 

It  seems  likely  that  uniaxial  compression  along  certain 
crystallographic  directions  will  activate  the  bcc  to  hep 
transformation.  The  present  search  for  possible  uniaxial  com- 
pression transformations  is  limited  to  compressions  of  the  tri- 
clinic cell  of  Fig.  5.2(b). 

This  cell  attains  hexagonal  symmetry  by  uniaxial  compres- 
sion in  the  ToOl]  direction;  then  the  70°30  angle  between 
sides  in  the  basal  plane  becomes  60°.  Due  to  the  symmetry  of 
cubic  crystalS^  six  different  planes  exist  which  contain  a tri- 
clinic cell  like  that  of  Fig.  5.2(b).  However,  there  are  only 
six  directions  of  uniaxial  compression  which  will  produce 
hexagonal  symmetry.  Planes  which  contain  triclinic  cells  and 
directions  of  uniaxial  compression  which  transform  the  cell  to 
hep  are:  • 


Plane 

Direction 

(lIO) 

[001]  ,[ooI] 

(110) 

[OOl]  ,[ooI] 

(011) 

[ioo],[Ioo’] 

(Oil) 

[loo]  ,[ioo  ] 

(101) 

[oio],[oIo  ] 

(lOl) 

[oio],roio  ] 

The  relative  displacement  required  to  obtain  the  60 
angle  is  a/Sp  = 0.816  = relative  volume  for  uniaxial  strain. 

To  produce  this  compression  in  iron  requires  a shock  stress  near 
400  kbar.  The  actual  transformation  occurs  at  130  kbar,  which 
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k. 


suggests  that  some  other  process  is  active  at  lower  stresses  or 
some  local  expansion  occurs  laterally. 

5.2.3.  Fraction  Transformed 

Theories  of  martensitic  transformations  have  progressed 
to  the  point  where  the  volume  fraction  transformed  can  be  ex- 
pressed as  a function  of  time  for  constant  temperature  and  pres- 
sure or  as  a time-independent  function  of  temperature  for 

47 

constant  pressure.  Transformations  of  the  first  kind  are 
called  "isothermal,"  the  latter  are  called  "athermal."  Iso- 
thermal data  show  that  reaction  rate  decreases  rapidly  during 

the  course  of  transformation.  This  has  been  explained  by 
46 

Fisher  as  a consequence  of  repeated  subdivision  of  grains  by 
succeeding  generations  of  martensitic  plates.  Plates  growing 
from  any  initiation  site  are  stopped  at  grain  boundaries  or 
other  plates  so  that,  as  a grain  subdivides,  the  transformed 
volume  growing  from  each  nucleation  site  diminishes  as  in 
Fig.  5.4.  If  nucleation  sites  are  activated  at  a constant  rate, 
and  if  this  rate  is  slow  compared  to  individual  growth  rate, 
then  it  necessarily  follows  that  volume  transformation  rate 
diminishes  with  time;  the  amount  of  material  transforming  de- 
pends on  the  amount  of  material  available  to  transform.  The 

"final"  state  reached  in  such  experiments  has  been  observed  to 

47 

contain  a residue  of  phase  1. 

In  constant  pressure  athermal  studies  a sample  is  heated 
above  the  martensitic  transition  temperature,  M^,  and  rapidly 
cooled  to  a reference  temperature  below  M,p.  It  is  then  examined 


J 


1 


1 

I 

: j 


FIGURE  5.4  SUBDIVISION  OF  A SINGLE  GRAIN  BY  MARTENSITIC  PLATES. 
A AND  B ARE  FIRST  AND  SECOND  GENERATION  PLATES, 
RESPECTIVELY. 
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(Of  (ho  omouiU  of  mar  tioniii  to  formod.  Tlic  amount  of  martensite 
formeii  has  been  found  proportional  to  the  difference  between 
transformation  temperature  anti  the  lower  reference  temperature 
ami  independent  of  time,  within  tlie  resolution  of  sucli  experi- 
ments. 

S . 2 . 4 . Fract ion  Trans  formed  in 
Athe  r ma 1 Trans  forma  t ion 

Athermal  martensite  appears  to  grow  from  a finite  number 

4 

of  nuc.leation  sites,  wl;ich  are  eventually  exliausted.  Magee 
has  been  successful  in  explaining  athermal  results  by  assuming 
the  number  of  riew  plates  in  the  unt  rans  formed  volume  to  be  pro- 
portional to  increments  in  "driving  force,"  G.,  ^ - G^, 

where  G.,  and  G^  arc  the  Gibbs  etuM:gie:'.  for  phases  1 and 
2,  respectively,  at  ttie  existing  pressure  and  temperature.  He 
has  proposed  that  the  number  of  new  sites  per  unit  volume,  dN, 
available  to  transform  is 

I 

dN  = -K(l-f')  d(G._,-G^)  (5.1) 

1 

wliere  K is  a positive  cons.tant  and  f is  volume  fraction  of 
mar  tens i te . 

Tills  equation  contains  the  implicit  assumption  that  tlie 

I 

final  state  is  tlie  state  for  which  f = I,  whereas  oliservat  ions 
show  that  f = I is  never  quite  reached. 

If  eacli  cluster  of  plates  contains  nearly  tlie  same  volume 
of  martensite,  t!ien  the  volume  fraction  transformed  is  propor- 
tional to  tlie  number  of  new  plates  formed  per  unit  volume: 
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df  = Vp  dN  (5.2) 

where  Vp  is  the  constant  average  plate  volume  of  a cluster. 

46 

This  treatment  of  plate  volume  is  different  from  Fisher's  and 
more  in  accord  with  athermal  experiments. 

An  expression  for  1 - f is  obtained  by  combining 
Eqs.  (5.1)  and  (5.2)  and  integrating,  which  gives 

1 - f'  = exptKVp(G2-Gj^-A)  J (5.3) 

where  A is  the  value  of  G2  - at  the  onset  of  transforma- 
tion. Magee  demonstrated  that  Eg.  (5.3)  or  some  equivalent 
accurately  fits  almost  all  results  from  athermal  martensitic 
rapid  cooling  studies. 

t 

The  volume  fraction  f is  equal  to  Pq/p2  times  the 
mass  fraction  f.  Since  Pg/p2  ^•‘5  nearly  1,  Eq.  (5.3)  is 
expected  to  be  valid  for  mass  fraction  data.  In  fact,  by  using 

I 

similar  assumptions,  Eq.  (5.3)  can  be  derived  with  f replaced 
by  f . 

The  search  for  a physical  model  to  reproduce  Eq.  (5.3) 
is  certainly  a major  problem  in  the  physics  and  metallurgy  of 
martensitic  transformations.  The  technological  importance  of 
these  transformations  in  commercial  materials  enhances  interest 
in  understanding  such  phenomena. 

5.3.  Nucleation  of  a New  Phase 

Various  types  of  nucleation  sites  are  discussed.  A new 
concept  is  described  in  which  nucleation  sites  are  "frozen-in" 
by  a sudden  step  in  pressure. 
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5,3.1.  Nucleation  Sites 

Energy  considerations  have  led  to  the  widely-accepted 
belief  that  nucleation  of  the  second  phase  occurs  on  defects 
such  as:  (1)  impurities,  (2)  grain  boundaries,  and  (3)  disloca- 
tions. Calculations  of  energy  to  nucleate  on  various  defects 
show  that  a dislocation  is  most  important. 

If  nucleation  occurs  on  defects,  it  is  reasonable  to 
assume  that  not  all  the  sites  are  activated  by  the  same  driving 
force.  An  increase  in  active  nucleation  sites  with  increasing 
driving  force  provides  a plausible  explanation  of  the  observed 
mixture  of  phases  for  final  states  in  athermal  martensitic  and 
certain  shock  experiments.  A satisfactory  theory  of  transforma- 
tion would  provide  a physical  mechanism  or  mechanisms  for  such 
an  increase,  but  such  has  not  yet  appeared. 

Twins  have  been  suggested  as  nucleation  sites,  but  no 
detailed  energy  calculations  have  been  made.  Since  twins  can  be 
formed  from  interacting  dislocations,  one  might  assume  that 
energy  requirements  for  nucleation  at  a twin  are  of  about  the 
same  magnitude  as  for  dislocations;  in  that  case,  twin  surfaces 
would  be  important. 

In  some  materials,  like  iron  and  steel,  twinning  is  the 

dominant  mechanism  for  plastic  yielding  under  rapid  compression. 

In  fact,  the  observed  rapid  increase  in  density  of  twins  as  the 

transformation  stress  in  iron  is  approached  has  led  German,  et 
30 

al.  to  suggest  twins  as  nucleation  sites  for  the  hep  phase. 

Consideration  of  the  bcc  lattice  shows  that  only  a 
small  change  is  required  to  obtain  hep  structure  along  the 
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plane  of  the  twin.  A crystallographic  model  illustrating  a 
twinned  bcc  lattice  is  shown  in  Fig.  5.5(a)  with  some  irregular 
hexagons  outlined.  The  lattice  is  viewed  as  stacked  (110) 
planes.  A small  shear  of  the  twin  plane  in  the  [ill]  direction 
results  in  a layer  of  symmetric  hexagonal  cells  along  the  twin 
plane  as  shown  in  Fig.  5.5(b).  A semicoherent  interface  exists 
between  the  two  phases.  The  shear  mechanism  is  similar  to  the 
one  given  in  Fig.  5.3.  It  appears  feasible  that  the  nucleation 
of  the  hep  phase  can  occur  on  twins. 

A possibility  which  has  not  been  previously  investigated 
is  that  equilibrium  embryos  of  the  second  phase,  which  always 
exist  as  a result  of  statistical  fluctuations,  may  be  "frozen- 
in"  by  sudden  application  of  sufficient  pressure  to  bring  the 
material  into  the  stability  field  of  the  second  phase.  This 
type  of  nucleation  would  apply  particularly  in  a shock  wave.  It 
is  described  in  detail  in  the  next  section. 

5.3.2.  Nucleation  Due  to  Rapid 
Application  of  Stress 

Consider  here  only  material  in  a lattice  that  is  perfect 
except  for  thermal  fluctuations.  In  an  equilibrium  state  of  the 
first  phase  there  is  a distribution  of  embryos  of  the  second 
phase  due  to  normal  statistical  fluctuations.  Figure  5.6(a) 
shows  the  number  of  embryos  in  stable  phase  1 as  a function  of 
the  number  of  atoms,  n,  in  each  embryo.  These  embryos  are 
created,  grow,  and  shrink  through  fluctuations.  For  this  case, 
where  phase  1 is  stable,  energy  to  grow  increases  monotonically 
as  shown  by  curve  A in  Fig.  5.6(b)  where  62^^  > 0.  If 


I 
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(a)  CROSS  SECTION  OF  A TWIN  (b)  SMALL  SHEAR  IN  THE  [TTi] 

IN  BCC  IRON  AS  VIEWED  ON  DIRECTION  IN  (112)  TWIN  PLANE 

THE  (110)  PLANE  RESULTS  IN  HEXAGONAL 


SYMMETRY  ALONG  TWIN  PLANE 


FIGURE  5.5  SHEAR  MECHANISM  IN  BODY-CENTERED-CUBIC  IRON  TO  OBTAIN 
HEXAGONAL  SYMMETRY  ALONG  A TWIN  PLANE. 
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(o)  DISTRIBUTION  OF  PHASE  2 
EMBRYO  SIZES  IN  STABLE 
PHASE  I 


AW 


AW, 

AW2 


FIGURE  5.6  SPHERICAL  EMBRYOS  IN  STABLE  PHASE  1.  CURVE  A IS  FOR 
THE  STABLE  PHASE  1 FIELD  WITH  G2  - G.,  = G21  > 0.  CURVE 
B IS  FOR  LOW  PRESSURE  SUFFICIENT  TO  MAKE  THE  PHASE  2 
FIELD  STABLE  WITH  G2^  (n/ ) < 0.  CURVE  C IS  FOR  HIGHER 

PRESSURE,  P2,  FOR  WHICH  Gj^  (n2*)<G2i  In,*)  <0. 
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sufficient  stress,  is  applied  to  make  the  field  for  phase 

2 stable,  G22  changes  sign  and  curve  B of  Fig.  5.6(b)  re- 

* 

suits.  For  this  case,  an  embryo  for  which  n < n^^  requires 

it 

energy  to  grow;  those  with  n^^  > n^^  will  continue  to  grow  be-  • 

cause  energy  of  the  system  is  diminished  by  growth.  Embryos  of 

_ ir 

size  n > n are  called  nuclei. 

Increasing  stress  to  P2  produces  energy  curve  C in 
Fig.  5.6(b).  In  this  case,  the  driving  force  G22^(P2)  is  more 
negative  than  ^21^^!^'  energy  required  to  form  a stable 

nuclei  of  phase  2 is  reduced  from  to  AW2.  For  curve  C 

, * * * 

j a stable  nuclei  of  phase  2 contains  n2  atoms;  n2  < n^^. 

If  the  distribution  of  embryo  sizes  given  in  Fig.  5.6(a) 

does  not  change  when  pressure  is  increased  form  zero  to  P2f 

then  the  number  of  embryos,  N2,  which  become  stable  includes 

* 

all  nuclei  for  which  n > n2.  This  condition  applies  if  the 
time  required  to  apply  the  driving  force  is  less  than  the  time 
required  to  redistribute  embryo  sizes  to  conform  to  the  new 
state.  This  should  certainly  be  true  for  large  embryos  in  a 
solid  where  reformation  times  are  large  when  a shock  wave  passes. 

A relationship  between  driving  force  and  number  of  nucle- 
ation  sites  based  on  these  concepts  can  be  derived  using  nuclea- 
tion  theory.  Assume  the  generally  accepted  conditions  of 
nucleation,  namely,  that  a number  of  small  clusters  of  atoms, 
called  embryos,  exist;  embryos  are  internally  uniform  and  have 
the  same  structure  and  properties  as  the  final  phase  in  bulk 
form.  These  assumptions  about  embryos  leave  shape  and  size  as 
the  only  variable  parameters.  A geometrical  shape  is  adopted  by 
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the  embryo  which  minimizes  the  energy  of  the  system  for  the  for- 
mation of  a nucleus.  For  simplicity,  assume  that  n atoms  in 
an  embryo  form  a sphere.  When  the  spherical  embryo  is  formed, 
energy  of  the  entire  assembly  increases  by  the  amount  AW 
defined  by  the  expression^^ 

AW  = n G23^  + a n^^^  (5.4) 


where  V2  is  volume  per  atom  of  the  second  phase  and  o is  sur- 
face energy  per  unit  area  of  the  interface  between  the  two 
phases.  The  surface  energy  term  in  Eg.  (5.4)  makes  it  necessary 
for  G22  to  attain  a certain  negative  value  before  transforma- 
tion can  begin. 

Consider  now  the  stability  of  embryos  according  to  their 

* 

size.  Embryos  containing  n atoms  maximize  AW,  where 


n = — - 


3G 


(5.5) 


21 


Having  reached  this  size,  they  continue  to  grow. 

In  macroscopic  assemblies,  fluctuations  will  lead  to 

local  transitory  phase  transformations.  These  fluctuations  are 

also  responsible  for  establishing  a distribution  of  embryos  of 

different  sizes  within  a stable  phase  which  has  the  same  atomic 

53 

arrangement  as  the  new  phase.  Frenkel  treated  each  embryo  as 
a molecule  of  a particular  kind,  independent  of  the  others,  and 
randomly  present  as  a dilute  solution  in  phase  1.  The  thermo- 
dynamic Gibbs  potential  of  such  a solution  is  given  by  the 
expression 
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G = 


NiGj^ 


I 

n 


N. 


n 


^2" 


(36tiV^) 


o n 


2/3 


(5.6) 


+ kT 


(/nA  / N \ 

N,  log  I — 1+  I N log!  —) 

\V  n \V 


where  N^j  is  the  number  of  embryos  containing  n atoms  and 
' - 

N ~ Z ‘3®‘^°tes  the  total  number  of  species,  those  of 


n 

different  size  being  treated  as  molecules  of  different  kinds. 
Size  distribution  of  clusters  was  determined  by  maximizing 
Eq.  (5.6)  to  obtain  the  expression. 


= N exp  , (5.7) 

3 

where  N = Z is  total  number  of  atoms  per  cm  , k is  the 

n 

gas  constant,  and  T is  temperature. 

53 

According  to  Frenkel,  the  kinetics  of  transformation 
in  solids  does  not  differ,  in  principle,  from  kinetics  of  the 
condensation  process  considered  above,  except  for  certain  spe- 
cial features  connected  with  the  shapes  of  crystals.  Therefore, 
we  assume  that  Eq.  (5.7)  is  valid  for  solid  systems. 

If  the  distribution  of  embryos  has  insufficient  time  to 
change  from  its  initial  form  when  a large  stress  is  rapidly 
applied,  a certain  number  of  embryos  are  suddenly  found  to  ex- 
ceed the  critical  or  stable  size  for  the  new,  thermodynamic 
state,  becoming  "trapped"  on  the  side  of  the  energy  curve  favor- 
ing growth.  The  number  of  embryos  that  begin  to  grow  under  these 
conditions  can  be  calculated  form  the  distribution  function  of 
stable  embryos  for  old  and  new  states.  The  number  of  embryos 

**12  interval  (n^^,  02)  is  given  by  the  expression, 
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"2  r^2 

Ni2  - I An  N„  dn  , 


(5.8) 


where  is  a large  number  and  An  is  -‘qual  to  1. 

Substitution  of  Eq.  (5.7)  into  Eq.  (5.8)  and  integration 

* 00  * 
over  the  interval  (n  ,n  ) gives  the  number  of  clusters  N 

that  exceed  the  critical  size  for  a sudden  change  in  driving 


force : 


* r o’) 

N = N r - AW^(n) 

J * exp 

•'n  L *^^0  J 


dn  , 


(5.9) 


where  AW°(n)  is  given  by  Eq.  (5.4)  for  the  initial  stable  state 

■k 

ahead  of  the  shock  front,  and  n is  given  by  Eg.  (5.5)  for  the 

state  behind  the  shock  front.  We  are  not  able  to  integrate 

Eq.  (5.9)  analytically.  However,  we  can  differentiate  Eq.  (5.9) 

54 

with  respect  to  G.,,  using  Leibnitz's  rule,  which  results  in 


the  expression. 


35^7 


. N expT-  dn- 

‘^°21  L ''’^0  J 3BJ7 

- N exp  [-  1 

L "To  Jasil 


3377 


(5.10) 


^^21 


Since  the  integrand  and  n are  independent  of  G21' 


J*  _ r AW°(n*)  1 dn* 

. — *“N  exp  — |.rp  jp  • 

’21  L ^^0  J °^21 


(5.11) 
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Differentiating  Eg.  (5.5)  with  respect  to  62^^  results 


in  the  expression. 


dn*  _ _* 


2 '^''2 


’21 


V2  dG2, 


(5.12) 


Substituting  Eg.  (5.12)  into  Eg.  (5.11)  results  in  the  expression, 

dV.. 


dN* 


dG 


= Nn* 


21 


’21 


V2  dG2i 


exp 


AW°(n*) 

kT 


(5.13) 


Eliminating  n*  and  AW°(n*)  by  Egs.  (5.4)  and  (5.5)  results 
in  the  final  expression. 


dN* 


-327Ta^V2N 


dG 


21 


3G 


21 


2 ^^2 


G21  ^2  ^^^21 


X exp 


-16770^2 


'21 


1 - 


2G22^(P=0,T=295”K) 


3G 


21 


, (5.14) 


which  depends  on  a,  dV2/dG22^,  and  62^^.  This  egua- 

tion  is  significant  because  .it  establishes  a relation  between 
number  of  nucleation  sites  and  driving  force,  G21,  in  the 
stable  field  of  phase  2. 

To  calculate  values  of  dN*/dG2j^f  from  Eg.  (5.14)  re- 

guires  values  for  surface  energy  o,  volume  V2,  N,  driving 

force  G21'  ^^2'^^^21'  Values  for  a found  in  the  litera- 

2 

ture  vary  from  20  ergs/cm  for  a coherent  twin  interface  to 
200  ergs/cm^  for  an  incoherent  inter face. Examination  of 
Eg.  (5.14)  reveals  that  the  influence  of  the  exponential  function 
is  overriding  and  if  any  nucleation  is  to  occur,  the  argument  of 
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the  exponential  must  be  less  than  about  100.  This  implies  that 

2 

the  value  for  o must  tiot  exceed  about  20  ergs/cm  . Values  of 

the  right  hand  side  of  Eg.  (5.14)  were  therefore  calculated  for 
/ 

2 

0 = 20  ergs/cm  . Values  for  the  driving  force  which  are 

given  in  multiples  of  A,  were  obtained  from  the  two-phase 
equation  of  state  for  iron  given  in  Appendix  A.  The  value  for 

— IP  T>T 

A = -8.33  x 10  ergs/atoms  is  G2\  P ~ 130  kbar, 

T = 338°K,  which  is  the  state  for  onset  of  transformation  for 

the  loading  process.  The  value  for  02^^  (P=0 ,T=295°K)  = 6.85 
-14 

X 10  ergs/atoms,  which  defines  the  initial  distribution 
function  shown  in  Fig.  5.5(a),  was  obtained  using  the  approxima- 
tion = (V^-Vj^)  (P-P*^)  - (S2-Sj^)  (T-332°K)  . (See  Appendix  A.) 

-21  2 

A constant  value  V2  = 1.045  x 10  cm  /atoms  for  iron  at 
TL  o 

P = 130  kbar,  T = 338  K was  used  since  V2  varies  only  from 
1.045  X 10“^^  cm^/atoms  to  1.021  x 10~^^  cm^/atoms  over  the 
stress  range  130  kbar  to  200  kbar.  The  number  of  iron  atoms, 

N,  in  1 cc  is  8.48  x 10^^.  The  value  for  (2/V2) (dV2/dG22) 

= 2K^/(V2-Vj^)  = -150  gm/Mbar  cm^,  where  is  isothermal  com- 
pressibility of  phase  2 and  the  approximation  dG2j^  = AVdP 
was  used.  (See  Appendix  A.)  Since  3/A  = 3.4  x 10^  gm/Mbar  cm^, 

(2/V2)  (dV2/dG22^)  can  be  ignored. 

Inspection  of  values  given  in  the  last  column  of 
Table  5.2  reveals  that  as  the  magnitude  of  G22^  increases  so 
does  the  magnitude  of  dN*/dG2j^.  Since  dN*/dG22  for  G2^  = A 
is  so  small  compared  to  values  at  higher  driving  force,  nuclea- 
tion  is  initiated  for  values  at  a driving  force  between  A and 
2A.  A smaller  value  for  o would  force  the  initiation  of 
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nucleation  nearer  A,  but  the  adjustment  appears  unwarranted 
considering  uncertainties  in  o.  Values  of  n*  become  unphysi- 
cally  small  for  driving  forces  exceeding  2A.  The  value  of  this 
model  and  suggested  ways  to  improve  it  are  discussed  in  the  next 
section. 

If  because  of  cold  work,  grain  growth,  or  other  energy 
concentrating  processes,  a number  of  permanent  spherical  embryos 
of  phase  2 exist  in  the  stable  field  of  phase  1,  then  the 
energy  required  to  reach  critical  size  is  less  than  for  pure 
homogeneous  nucleation.  In  fact,  suppose  that  these  pre-existing 
sites  ate  of  size  nQ  and  that  normal  statistical  fluctuations 
created  a distribution  of  sizes  as  in  the  homogeneous  case;  then 
Eqs.  (5.4)  and  (5.5)  are  unchanged,  except  that  n + n^  is  sub- 
stituted for  n.  The  important  change  comes  in  that  N becomes 
the  number  of  pre-existing  sites  Nq  of  size  nQ.  Substitution 
of  Nq  for  N in  Eq.  (5.14)  results  in  an  expression  for 
dN*/dG2j^  for  the  case  fo  pre-existing  sites. 


5.4.  The  Quasistatic  Hugoniot 


Shock  data  described  in  Section  5.1  show  that: 


1.  Initial  transformation  rate  for  close-in  experiments 
exceeds  2 x lO^/sec  for  mass  fraction  of  phase  2. 

2.  Any  transformation  taking  place  behind  the  plastic 
II  shock  is  slow  since  there  is  no  evidence  of 
changes  in  stress  behind  the  shock  front. 

3.  Hugoniot  P-V  states  for  130  to  200  kbar  are 
vastly  different  from  those  predicted  by  equilib- 
rium hydrostatic  thermodynamics  and  measured  values 
of  dP/dT. 
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Discussion  of  points  1 and  2 is  deferred  to  Chapter  6.  Concern  | 

here  is  with  point  3 and  the  quasistatic  Hugoniot  data  in  the  | 

t] 

mixed  phase  region  shown  in  Fig.  4.7.  ; 

It  can  be  seen  from  Fig.  4.7  that  changes  in  V are  !' 

nearly  linearly  related  to  changes  in  pressure  for  P between 
130  and  200  kbar.  If  these  measured  points  deviate  from  equi- 
librium because  transformation  stops  before  completion,  then  the 
mass  fraction,  f(P,T),  can  be  calculated  for  each  point  on  the 
Hugoniot  above  the  transition  point  according  to  the  lever  rule, 

V (P,T)  - V(P,T) 

1 - f = , (5.15)  ; 

V2(P,T)  - Vj^(P,T) 

! 

where  f is  the  mass  fraction  of  phase  2,  V(P,T)  is  specific  I 

volume  of  the  mixture,  and  Vj^(P,T)  and  V2(P,T)  are  specific 
volumes  of  phases  1 and  2,  respectively.  As  V approaches 
^2  value,  small  systematic  uncertainties  in  ^2' 

and  V lead  to  large  uncertainties  in  1 - f:  i 

6(;-f)  _ '*''2  6V  

1-f  ~ ^2~^1  * 

Values  of  f calculated  from  Eq.  (5.15)  and  the  meas- 
ured Hugoniot  points  shown  in  Fig.  4.7  are  represented  in  ' 

Fig.  5.7,  where  ln(0.93-f)  is  plotted  versus  P-P  . This  odd 
ordinate  was  chosen  because  it  was  found  that  subtracting  0.07 
from  1-f  produced  a reasonable  fit  of  the  data  to  a straight 
line,  as  shown.  This  deviation  of  1 - f from  zero  at  the 
asymptote  may  be  due  to  error  in  estimating  the  position  of  the 
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phase  2 boundary  or  to  incomplete  transformation  for  large 
TL 

values  of  P - P . The  data  shown  are  insufficient  to  distin- 
guish between  these  two  possibilities,  each  of  which  is  repre- 
sented in  Eq.  (5.16).  The  third  term  on  the  right  hand  side  of 
that  equation  is  negligible. 

The  second  phase  surface  is  calculated  from  a two-phase 
equation  of  state  which  is  based  on  data  of  Mao,  et  al. (See 
Appendix  A.)  They  report  an  uncertainty  in  initial  volume  for 


the  second  phase,  Vq2,  of  .0011  cm^/gm.  This  uncertainty 

was  reduced  approximately  50  percent  using  a value  of  Vq2 

consistent  with  x-ray  measurements  of  V2-Vj^  made  at  stresses 

near  130  kbar  and  reported  by  Mao,  £t  al . The  difference, 

V2-V,  obtained  from  this  equation  of  state  and  measurements  by 

Barker  and  Hollenbach^^  goes  from  0.00056  to  0.0004  cm^/gm 

for  stresses  from  204  to  304  kbar,  which  suggests  that  the 

second  phase  surface  and  the  Barker  and  Hollenbach  data  agree 

within  uncertainties  of  the  experiments  and  accuracy  of  the 

equation  of  state  for  the  second  phase. 

TL 

since  P - P is  nearly  proportional  to  62^^  for  iron 
(see  Appendix  A) , an  equally  good  fit  is  obtained  by  plotting 
ln(.93-f)  versus  G2j^«  The  equation  of  the  line  so  obtained  is 

0.93  - f = exp[0  (Gt^^-A)]  , (5.17) 

where  9 = 6,444  gm/Mbar  cm^  is  determined  by  least  squares, 

-5  3 

and  A = 8.7  X 10  Mbar  cm  /gm  is  G2j^  at  the  transformation 
state  (p'^^,t'^^)  . The  differential  of  Eq.  (5.17)  is 
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df  = -e  (.93-f)  d(G2j^)  , (5.18) 

47 

which  corresponds  closely  with  the  relation  found  by  Magee  to 

describe  athermal  results  for  the  gamma  to  alpha  transformation. 

Eg.  (5.3).  Although  there  is  no  obvious  connection,  Magee  found, 

for  carbon  alloys  with  less  than  1.2  percent  carbon,  a value 

for  9 which  is  14  times  larger  than  0 from  Fig.  5.7. 

Due  to  the  large  uncertainties  in  1 - f for  large  values 

of  P shown  in  Fig.  5.7,  one  might  consider  data  only  for 

stresses  less  than  or  equal  to  204  kbar.  These  data  are  well 

fitted  by  a linear  relation  between  S,n(l-f)  and  G21'  They 

yield  a value  of  0 equal  to  4,048  gm/Mbar  cm^,  less  than 

two-thirds  the  previous  value. 

Equation  (5.17)  can  be  interpreted  in  the  following  way; 

G2]L  driving  force  which  causes  the  transformation  to  proceed. 

It  must  exceed  a threshold  value.  A,  before  the  transformation 

is  initiated.  For  unknown  reasons  the  transformation  proceeds 

only  as  G2j^  is  increased.  If,  for  example,  nucleation  sites 

exist  which  are  activated  at  different  stress  levels,  increases 

in  P produce  increases  in  622^,  more  nucleation  sites  are 

activated,  and  the  transformation  proceeds  incrementally.  The 

calculation  of  nucleation  sites  described  in  Section  5.3  provides 

a detailed  model  of  such  a situation,  Eq.  (5.14).  In  that  case, 

however,  the  proportionality  parameter  0 is  not  constant  but 

varies  over  a wide  range  of  values  in  the  mixed  phase  region. 

If  this  transformation  is  martensitic  and  transformation  occurs 

—8  3 

with  constant  average  plate  volume,  Vp,  of  10  cm  , then 
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'I 

I 

i 

I 


values  of  0 are  obtained  by  multiplying  Vp  by  entries  in  the 

-2  3 

last  column  of  Table  5.2.  This  gives  6 = 7 x 10  /Mbar  cm 

14  3 

for  G2j^  = A and  0 = 1.2  x 10  gm/Mbar  cm  for  G22  = 2A. 

These  values  are  fat  from  observed  values,  which  is  not  surpris- 
ing considering  the  unreality  of  the  basic  assumption  that 
nucleation  is  occurring  in  the  homogeneous  lattice. 

Values  for  0 calculated  from  Eg.  (5.14)  come  much 
closer  to  the  measured  value  when  the  homogeneous  model  is  modi- 
fied by  assuming  spherical  pre-existing  embryos  of  phase  2.  For 

—8  3 

this  case,  if  martensitic  plate  volume  Vp  = 10  cm  is  1 per- 
cent of  the  grain  volume  for  0. 1-mm-diameter  grains,  Nq  can 
be  calculated  by  dividing  values  of  dN*/dG2j^  in  the  last  column 
of  Table  5.2  into  -N/Vp  0 = 6.4  x 10^^.  This  gives  Ng  = 8.3 
X 10^®/cm^  for  G22^  = A,  Ng  = 5.3  x 10^^/cm^  for  G22  = 2A, 

Q O 

and  Ng  = 1.28  x 10  /cm-"  for  ^21  =3A.  These  values  for  Ng 
approach  seemingly  realistic  values  for  pre-existing  sites  since 
the  number  of  twins  required  to  account  for  all  the  plastic 

strain  in  shocked  iron  at  130  kbar  was  inferred  from  Johnson 

56  7 3 

and  Rhode  to  be  about  10  /cm  . 

The  above  calculations  show  how  a relation  between  dN 
and  dG22  be  established,  and  although  a detailed  model  to 
explain  the  form  of  Eq.  (5.17)  has  not  been  produced,  its  simi- 
larity to  the  Magee  equation  strengthens  the  link  between  the 
alpha  to  epsilon,  shock-induced  transformation  and  the  athermal, 
martensitic,  gamma  to  alpha  transformations  in  iron.  The  calcu- 
lation also  suggests  a basis  for  understanding  of  both  athermal 
gamma  to  alpha  and  shock-induced  alpha  to  epsilon  transformations. 
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The  basic  concept  of  freezing-in  statistical  fluctuations  appears 
sound  and  avoids  the  theoretical  problems  of  nucleation,  wherein 
growth  of  embryos  is  too  slow  to  enable  the  transformation  to 
proceed  as  observed.  Since  the  weight  of  observations  on  ather- 
mal  martensite  suggests  that  nucleation  occurs  at  lattice  imper- 
fections, this  calculation  should  be  repeated  for  the  much  more 
difficult  problems  of  embyos  formed  about  such  sites.  The 
simple  modification  to  homogeneous  nucleation  described  above 
suggests  that  such  calculations  will  be  in  the  right  direction 
to  bring  about  agreement  between  measurements  and  theory.  Con- 
siderations are  also  required  on  the  effects  on  the  transforma- 
tion of  (1)  strain  in  the  lattice  to  accommodate  the  nuclei, 

(2)  use  of  the  hydrostatic  Gibbs  function  in  problems  dealing 
with  solids,  and  (3)  use  of  bulk  values  of  surface  energy  on 
small  surfaces  with  large  curvature.  If  these  modifications  can 
be  accomplished,  it  may  well  turn  out  that  the  too  rapid  varia- 
tion of  dN/dG22  with  ^2\'  shown  in  Table  5.2,  vanishes  and 
that  a reasonable  theoretical  basis  for  the  observations  is  pro- 
duced. 
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CHAPTER  6 

KINETICS  OF  TRANSFORMATION 
20 

Hone  and  Duvall  treated  wave  propagation  in  a phase- 
transforming  material  by  assuming  that: 

1.  Phases  do  not  separate  mechanically;  i.e.,  particle 
velocity  is  common  to  both  phases. 

2.  Pressure  and  temperature  are  common  to  both  phases. 

3.  Interfacial  energy  is  a small  fraction  of  the  total 
energy. 

4.  Mass  transfer  is  irreversible  and  df/dt  is  a 
known  function  of  the  state  variables. 

These  assumptions  can  be  incorporated  in  the  flow  equations  to 

provide  a theoretical  basis  for  calculating  evolution  of  a shock 

wave  in  the  transforming  material,  including  the  effects  of 

transformation  kinetics.  Their  principal  result  can  be  expressed 

in  the  form: 


dP 

dt  ■ 

dp  df 

^11  dt  ^12  3t  ' 

(6.1) 

dT 

Jt  = 

dp  df 

®21  at  ^ ®22  eft  ' 

(6.2) 

df 

3t  = 

’i'(G2i,e)  , 

(6.3) 

where  a^^j^  is  frozen  sound  speed  at  mass  fraction  f;  the  other 
a^^j  are  functions  of  V,T,P,f.  They  assumed  for  calculational 
purposes  a simple  relaxation  function  for  'V{G2iii)i 
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(6.4) 


wn  M'e  i;;  t ; elaxation  time. 


27  29 

Andrews  ' modified  t to  a form  more  obviously  re- 


lated to  the  Gibbs  function: 


= VJI  = 


|a|t. 


, 0 < f®^  < 1 , 


(6.5) 


where  J is  the  Jacobian,  J = 3 (V,E) /3  (P,T)  ^ = “Cpj^VKg^,  and 
jA|  is  the  matrix  of  coefficients  of  differential  equations 
describing  the  mixed  phase  region.  (See  Appendix  D.) 

Relaxation  time,  t,  is  not  uniquely  defined  and  must 

20 

be  associated  with  a specific  model  for  f®*^.  Horie  and  Duvall^ 
assumed  f®^  to  be  defined  by  states  on  the  equilibrium  PVT 
surface  in  the  mixed  phase  region  at  the  same  volume  and  tempera- 
ture as  the  considered  state.  Figure  6.1  shows  a cross  section 
of  the  PVT  sui face  at  a constant  T and  two  possible  transient 
states  of  the  transforming  mixture,  A and  B.  State  A pro- 
jects to  a volume  V^  on  the  isotherm,  giving 


T 

v;  - v„ 

.fSq  _ _1 A 

T T 

- V 

^2  '^1 


(6.6) 


which  is  a function  of  pressure  and  temperature.  State  B is  a 
more  completely  transformed  state  at  the  same  temperature.  For 
this  case  Vg  < V^  and  f®^  = 1. 

Andrews  assumed  f®*^  to  be  defined  by  states  on  the 
equilibrium  PVE  surface  for  the  same  volume  and  energy  as  the 
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There  are  three  types  of  existing  experimental  data 
which  measure  different  aspects  of  shock  evolution  in  iron  and 
allow  evaluation  of  relaxation  time:  (1)  amplitude  of  stress 
behind  the  plastic  I shock  as  a function  of  sample  thickness, 

(2)  rise  time  of  the  plastic  II  shock,  and  (3)  residual  metal- 
lurgical effects  when  the  plastic  II  shock  is  rapidly  diminished 
by  relief  waves.  Values  of  and  bounds  for  x for  these  differ- 
ent types  of  data  are  discussed  separately  in  the  next  three 
subsections  and  are  shown  to  agree  within  about  50  percent. 

6.1.1.  Decay  of  Stress  Associated  ' 

with  the  Plastic  I Shock 

The  plastic  I shock  decays  in  amplitude  as  it  propagates. 

If  material  strength  is  ignored,  plastic  I amplitude  is  expected 
to  be  near  final  driving  amplitude  at  positions  near  the  impact 
boundary.  By  making  the  essential  assumption  that  phase  1 
material  is  shocked  to  final  driving  stress,  at  the  impact 

surface  and  that  the  plastic  I shock  velocity  is  constant, 

Eqs.  (6.1)  and  (6.4)  can  be  used  to  derive  an  approximate  dif- 
ferential equation  for  rate  of  decay  of  the  plastic  I shock 
amplitude : 
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dP 

ar 


df 


2V 


2 dt 


(6.8) 


where  AV  = V2-Vj^. 

Setting  f = 0 in  Eq.  (6.4)  and  assuming  constant  tem- 
perature and  linear  P-V  relations  for  phase  1 and  phase  2 
near  the  mixed  phase  boundaries,  they  obtained  the  relation. 


P = P*^^  + (P^-P*^^)  exp 


,TL 


X 


0 ^ < 1 , 


(6.9) 


where  P is  transition  stress,  P is  driving  stress,  x is 
sample  thickness,  and  U2  is  constant  plastic  I shock  velocity. 
Equation  (6.9)  was  found  compatible  with  data  in  Fig.  4.3,  assum- 
ing = 0.05  ysec  to  be  constant  for  0 < x < 1 mm  and  a 
final  driving  stress  of  201  kbar.  This  value  represents  an 
approximate  upper  bound  for  since  a 20  percent  increase 
in  its  value  is  incompatible  with  the  data,  while  effects  of  a 
decrease  are  undetectable. 

Barker  and  Hollenbach^^  found  that  = 0.17  ysec  was 
required  in  Eq.  (6.9)  to  explain  their  data  on  plastic  I stress 
decay  for  equal  propagation  distances  but  different  final  driv- 
ing stresses. 

Substituting  Eq.  (6.5)  into  Eq.  (6.8)  results  in 


P = p"^^  + (P^-p"^^)  exp 


1 


a^^AV^Jx 


0 < f®'^  < 1 , 


(6.10) 


27  2Q 

winch  describes  plastic  I decay  according  to  Andrews'  model.  ' 

2 

For  iron,  the  term  -CpAV  in  |A|  (see  Appendix  D)  exceeds  the 
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others  by  an  order  of  magnitude.  This  gives  J/|a|  = 
V^Kg/(<iV)^  = V^/(a^j^AV^)  which  makes  Eg,  (6.10)  identical  to 


Eg.  (6.9).  This  shows  that,  under  the  assumptions  used  to 
obtain  Eg.  (6.9),  = i2  Cor  iron.  This  identity  is  not  a 

general  result,  and  in  transformations  where  terms  other  than 

2 

-CpAV  dominate,  the  value  of  jAj  will  produce  different 
values  for  and  t^. 

E'ju^tions  (6.9)  and  (6.10)  are  strongly  dependent  on  the 


basic  assumption  that  the  shocked  phase  1 material  remains  in 
phase  1 and,  at  the  impact  surface,  reaches  the  driving  stress 
at  a point  on  the  metastable  or  extended  phase  1 surface. 

This  assumption  may  be  invalidated  by  the  inability  to  prepare 
smooth  microscopic  surfaces.  Even  the  best  finely-lapped  and 
polished  surfaces  contain  microvoids  which  require  closing 
before  stress  at  the  impact  surface  can  be  sustained.  If  the 
effective  driving  stress  were  to  be  thus  reduced,  T2  in  the 
above  equation  could  be  increased  without  violating  the  data. 

Data  of  Fig.  4.3  show  plastic  I first  decreasing  as  x 
increases,  then  increasing,  then  decreasing  again.  This  behavior 
might  arise  from  inaccuracies  in  measurements  which  have  not 
been  fully  accounted  for,  or  from  other  effects  such  as  shear 
strength  associated  with  precursor  decay  or  relaxation  in  the 
plastic  I wave  and  behind  it. 


Rise  times  of  0.2-0. 3 jisec  for  plastic  II  shocks  have 
23 

been  reported  elsewhere;  these  are  consistent  with  the  present 
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measurements,  which  gives  a value  of  0.18  + 0.02  ysec  for  the 
rise  time.  These  rise  times  are  used  here  to  evaluate  relaxation 
times  from  calculated  steady  shock  wave  profiles. 

Swan,  et  al . have  shown  that  transient  effects  in  shock 
fronts  ultimately  decay,  leaving  a steady  shock  profile.  If  it 
is  assumed  that  measured  profiles  of  the  plastic  II  shock  front 
are  steady,  and  that  the  locus  of  states  within  these  shocks 
lies  along  the  appropriate  Rayleigh  line,  it  follows  that  an 
estimate  of  relaxation  time  can  be  obtained  by  comparing  measured 
profiles  with  calculated  steady  profiles,  assuming  that  trans- 
formation kinetics  dominates  the  shock  transition  process. 

57  58 

Techniques  for  calculating  profiles  are  described  elsewhere.  ' 
The  technique  used  here  is  described  in  Appendix  D.  Values  of 
relaxation  time  inferred  in  this  way  should  be  greater  than  true 
values  because  effects  of  viscosity  act  to  increase  rise  time. 

Figure  6.2  shows  calculated  temperature-independent  pro- 
files for  = 0.05  Msec  and  i = 0.1  Msec.  The  relaxation 
time,  I = 0.05  nsec,  in  the  rate  equation  of  Eq.  (6.4)  was 
required  to  obtain  a wave  profile  with  rise  time  near  0.2  Msec. 
Width  of  the  shock  wave  as  used  here  is  defined  as  rise  time 
times  laboratory  velocity,  ^2  ~ *^1  " 3.36  mm/Msec.  A practical 
definition  of  5 percent  to  95  percent  of  maximum  stress 
amplitude  was  used  for  rise  time  and  shock  width  determinations. 

6.1.3.  Metallurgical  Data 

Residual  hardness  for  iron  which  has  been  shocked  into 
the  epsilon  phase  by  successive  plastic  I and  plastic  II  waves 
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is  several  times  greater  than  for  that  which  ha.s  ntit  been  tran.-. - 
formed.  Thickness  of  the  plastic  II  shock  front  has  been  esti- 
mated by  some  workers  from  post  shock  hardness  measurements  in 

which  the  plastic  II  wave  is  quenched  by  reflection  of  elastic 

24  25 

precursor  and  plastic  I wave  from  a free  surface.  ' 

Shock  width  can  be  related  to  the  hardness  transition 

[ zone  thickness  measured  in  such  experiments  in  the  following 

! 

way:  Assume  the  interaction  between  the  oncoming  plastic  II 

wave  and  reflected  plastic  I wave  to  be  as  shown  in  Fig.  6.3. 

Then 

h2~hT  , , , , 

r-^  = {hj^-h2)  (I/R2  - l/R^)  + (h2-h5)(l/U2  + I/R2)  r (6.11) 

^2 

where  velocities  are  Lagrangian,  shock  width, 

hj^-h2  is  distance  from  the  free  surface,  h2-hj  is  width  of 

t I 

the  transition  band,  U2  is  plastic  II  wave  velocity,  R2  is 
velocity  associated  with  the  relief  wave  that  brings  the  stress 

I 

behind  the  plastic  II  shock  to  130  kbar,  and  is  the  lead- 
ing elastic  relief  wave  velocity.  In  deriving  this  expression 
the  following  assumptions  were  made:  (1)  wave  velocities  are 
constant  through  the  region  of  wave  interactions;  (2)  the  plas- 
tic II  shock  is  steady  and  can  be  represented  as  two  parallel 
lines  in  the  h-t  diagram;  and  (3)  the  first  term  on  the  right 
hand  side  of  Eg.  (6.11)  is  smaller  than  the  second  term.  Slopes 
of  the  various  lines  are  indicated  in  the  h-t  diagram  as  the 
reciprocal  of  their  wave  speeds. 
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Average  values  of  = 2.53  mm  and  h2~hg  = 0.47  mm 

were  found  for  a 190-kbar  shock  in  12-mm-thick  samples  of 

25  ' ' 

iron.  If  U2  = 3.578  mm/pseo,  R2  = 5.69  mm/|iSf?c,  and 

I 

= 6.58  mm/psec,  then  h2-h2  = .98  mm  which  is  equivalent 
to  a rise  time  of  0.27  ysec  for  the  plastic  II  wave  front. 

This  value  is  within  the  range  of  observed  values  of  rise  time 
described  in  the  preceding  section,  but  near  the  high  side. 


6.2.  Slow  Decay  of  the  Stress 
Behind  the  Plastic  I Shock 

TL 

It  was  noted  in  Chapter  4 that  P diminishes  slowly 
with  propagation  distance;  but  that  if  stress  jump  across  the 
plastic  I front  is  considered,  this  slow  decay  disappears. 

TL 

Therefore,  one  can  reasonably  infer  that  the  slow  decay  of  P 
is  due  to  precursor  decay.  The  situation  can  be  complicated  by 
wave  interactions  so  the  inference  is  not  conclusive. 

The  situation  can  be  clarified  by  describing  possible 
bounds  of  stress-particle  velocity  states  at  the  impact  boundary 
when  an  iron  sample  is  impacted  by  an  aluminum  projectile. 

Figure  6.4  illustrates  the  preusure-par t icle  velocity  plane; 
dashed  lines  represent  metastable  extensions  of  lower  pressure 
states  and  the  solid  lines  represent  equlibrium  Hugoniots.  The 
aluminum  cross  curve  represents  possible  states  at  the  impact 
boundary.  Point  A represents  the  maximum  attainable  stress  at 
the  instant  of  impact,  and  C represents  the  equilibrium  stress 
obtained  when  all  time  effects  have  disappeared.  The  problem  of 
kinetics  at  the  impact  surface  is  to  describe  how  fast  the  stress 
gets  from  A to  C.  The  time  required  for  stress  to  decay  from 
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FIGURE  6.4  STRESS-PARTICLE  VELOCITY  STATES  AT  THE  IMPACT  BOUNDARY  WHEN  AN 
IRON  SAMPLE  IS  IMPACTED  BY  AN  ALUMINUM  PLATE. 
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A to  B depends  on  the  elastic-plastic  transformation  rate, 
while  the  time  for  stress  to  decay  from  B to  C depends  on 
the  polymorphic  transformation  rate.  To  understand  the  rate  of 
decay  from  A to  C requires  an  understanding  of  both  transfor- 
mation processes.  However,  if  the  yielding  rate  is  much  faster 
than  the  polymorphic  transformation  rate,  the  polymorphic  trans- 
formation rate  is  the  limiting  and  important  one.  In  that  case, 
details  of  the  yielding  process  are  not  necessary  to  solve  most 
problems  of  interest. 

Both  precursor  and  plastic  I decay  data  show  that  both 
transformations  occur  very  rapidly.  However,  the  data  of  Chap- 
ter 4 are  not  sufficient  to  resolve  effects  of  the  kinetics  of 
yielding  on  the  evolution  of  shocks  in  iron.  Further,  experi- 
mental and  theoretical  work  in  this  area  are  required. 

6.3.  Discussion 

There  are  some  disparities  in  relaxation  times  inferred 
from  various  types  of  data  which  suggest  that  relaxation  time  is 
not  constant.  However,  the  relaxation  times  are  all  of  the  order 
of  0.1  ysec.  A lower  bound  of  initial  transformation  rate  in 
iron  shocked  to  200  kbar  of  10^/sec  can  be  inferred  from  this 
value  of  T and  Sq.  (6.4). 

The  theoretical  treatment  appears  adequate  for  data 
being  obtained.  Future  work  may  require  closer  examination  of 
its  basic  assumptions. 
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CHAPTER  7 

SUMMARY  AND  CONCLUSIONS 

Little  is  known  about  atomic  mechanisms  which  cause 
ultrafast  phase  transitions.  This  void  in  knowledge  is  partly 
due  to  lack  of  experimental  data.  Need  for  data  and  understand- 
ing of  the  transformation  process  is  the  reason  for  this  study. 

Iron  was  chosen  as  the  material  to  study  primarily  because  of  | 

the  wealth  of  existing  thermodynamic  and  equation  of  state  data.  i 

! 

Experimental  intent  of  this  study  was  to  measure  evolu-  i 

tion  of  the  plastic  I shock  in  polycrystalline  Armco  iron  when 
final  driving  stress  is  near  200  kbar.  The  most  significant 
experimental  result  is  that  little  or  no  variation  of  plastic  I 
wave  amplitude  occurs  for  propagation  distances  between  0.9  t 

and  6.35  mm.  This  implies  a relaxation  time  of  about  0.05  yseo 
or  less  for  initial  stages  of  the  alpha  to  epsilon  transformation. 

When  the  single  measurement  at  25.4-mm  propagation  distance  is 
included  with  close-in  data,  one  infers  a slow  variation  of 
transformation  stress  with  distance  superimposed  on  the  rapid 
decay  for  propagation  distance  of  less  than  1 mm.  Inference  of 
initial  decay  rate  depends  on  the  assumption  that  initial  com- 
pression at  the  impact  surface  is  entirely  in  the  alpha  phase. 

Further  results,  which  are  essentially  corrobative,  are 
that:  i 


J 
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1.  Elastic  precursor  amplitude  increases  as  sample 
thickness  is  decreased. 

2.  Transformation  stress  measured  in  a 25.4-mm-thick 
sample  is  131.4  + 3.3  kbar. 

3.  Relative  volume  behind  the  plastic  II  wave  is 
Vj/Vq  = 0.871  + .008  for  a stress  of  201  + 8.4  kbar. 

4.  Rise  time  in  the  plastic  II  shock  front  is 

0.18  + .02  usee. 

These  results  are  in  essential  agreement  with  those  of  earlier 
experiments  reported  in  References  36,  2,  2,  and  23,  respectively. 

In  addition  to  this  study  other  data  exist  which  relate 
to  kinetics  of  transformation.  They  are: 

1.  Rise  times  of  0.2-0. 3 ysec  for  the, plastic  II 
shock  measured  by  Novikov,  et  al . and  more  re- 
cently by  Barker  15  and  this  author . 

2.  Metallurgical  measurements  of  a hardness  transfor- 
mation zone  thickness  in  shocked  iron  by  Smith24 
and  more  recently  by  Smith  and  Fowler. 25 

3.  Slow  decay  of  the  plastic  I shock  in  iron  first 
observed  by  Minshall22  and  confirmed  by  this  work. 

Considering  all  but  the  slow  decay  of  plastic  I stress,  the  in- 
ferred relaxation  time  from  the  data  was  approximately  equal  to 
0.1  + .05  ysec.  This  gives  a lower  bound  for  the  initial  trans- 

7 

formation  rate  of  10  /sec  for  the  second  phase. 

The  slow  decay  in  the  plastic  I shock  amplitude  may  be 
due  to  decay  in  the  elastic  precursor  amplitude  for  iron  samples 
thicker  than  1 mm.  This  is  suggested  by  the  apparently  con- 
stant stress  jump  in  the  plastic  I shock.  Because  of  effects  of 
wave  interactions,  this  conclusion  is  not  firm;  experiments  with 
samples  treated  to  modify  precursor  structure  may  improve  under- 
standing of  this  phenomenon. 
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Based  on  the  assumption  that  Hugoniot  data  between  130 
and  250  kbar  are  for  partially  transformed  material,  the  lever 
rule  was  used  to  calculate  the  mass  fraction,  f,  of  phase  2 
which  is  experimentally  related  to  excess  stress  above  the  trans- 
formation stress.  These  data  show  that  f varies  exponentially 
with  G2j^,  the  difference  between  Gibbs  energies  of  the  bulk 
phases,  giving  the  expression 

1 - f = exp[9 (G2^-A) J , 

where  6 is  constant  and  A is  the  value  of  G2]^  at  onset  of 
transformation. 

Previous  workers  have  shown  that  such  a relationship 
describes  athermal  martensitic  transformations;  this  similarity 
strengthens  the  link  between  martensitic  transformation  and  the 
shock-induced  alpha  to  epsilon  transformation.  A possibility 
which  has  not  been  previously  investigated  is  that  equilibrium 
embryos  of  the  second  phase,  which  always  exist  as  a result  of 
statistical  fluctuations,  may  be  "frozen-in"  by  sudden  application 
of  sufficient  pressure  to  bring  the  material  into  the  stability 
field  of  the  second  phase.  A relation  between  number  of  "frozen- 
in"  nucleation  sites  and  driving  force  can  be  established;  this 
fact  suggests  a basis  for  understanding  both  athermal  transfor- 
mations and  the  shock-induced  alpha  to  epsilon  transformation. 

Better  theoretical  calculations  on  "frozen-in"  nuclea- 
tion sites  are  needed  which  include  strain  effects  and  surface 
energy  of  small  clusters  of  atoms;  errors  which  result  from  use 
of  the  hydrostatic  Gibbs  energy  for  solids  need  to  be  evaluated. 
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Experiments  to  date  have  not  revealed  the  operating 
mechanism  for  the  ultrafast  transformation  in  shocked  solids. 

If  the  alpha  to  epsilon  transformation  is  martensitic  as  sug- 
gested, then  experiments  with  different  grain  sizes  should  give 
different  kinetic  results  at  close-in  distances.  Single  crystal 
studies  would  be  useful  since  it  may  be  a shear-induced  trans- 
formation with  preferred  directions  and  planes.  If  the  trans- 
fv-rmation  is  nucleated  on  twins,  cold  rolling  to  eliminate  twin 
formation  prior  to  shock  loading  might  change  the  kinetics  of 
transformation.  Double  shock  experiments  in  which  first  shock 
amplitude  is  varied  in  order  to  vary  the  shock-induced  twin 
density  should  be  interesting. 

Surface  roughness  may  prevent  stress  on  the  metastable 
or  extended  phase  1 surface  from  reaching  the  driving  stress 
which  would  violate  a basic  assumption  in  Horie  and  Duvall's 
theory  of  plastic  I decay  ani  estimates  of  relaxation  time 
inferred  from  close-in  measurements.  A possible  experiment  to 
eliminate  surface  effects  would  be  to  deposit  within  a few 
nanoseconds,  at  some  depth,  enough  high  energy  electrons  to 
create  a shock  exceeding  the  transformation  pressure.  This 
would  not  be  a simple  experiment,  but  it  would  avoid  the  surface 
problem  since  the  shock  would  be  formed  internally. 
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APPENDIX  A 

THERMODYNAMICALLY  CONSISTENT  EQUATION  OF  STATE 

The  present  interest  in  an  equation  of  state  is  restricted 
to  describing  isotherms  and  Hugoniots  of  alpha  and  epsilon  iron. 
These  thermodynamic  paths  involve  very  small  portions  of  the 
equilibrium  surface  of  each  phase.  In  fact,  when  iron  is  shocked 
to  a stress  of  200  kbar,  the  total  temperature  rise  is  less 
than  100°K.  This  temperature  rise  will  displace  the  Hugoniot 
stress  by  less  than  3 percent  from  the  isothermal  stress  at 
the  same  volume. 

By  assuming  that  the  specific  heat  is  constant,  simple 

analytical  expressions  result  which  describe  the  equilibrium 

alpha  and  epsilon  surfaces.  Accurate  isothermal  compression 

data  have  been  reported  for  alpha  and  epsilon  iron  in  terms  of 

33 

the  Mucnaghan  equation.  The  equation  of  state  developed  below 

reduces  exactly  to  the  Murnaghan  equation  for  isotherms.  The 

accuracy  of  the  calculated  Hugoniot,  which  varies  only  slightly 

form  the  isotherm,  is  limited  principally  by  the  experimental 

error  of  + 2 percent  in  stress  for  the  measured  isotherms. 

The  choice  of  the  present  equation  of  state  rather  than 
27  28 

that  of  Andrews  ' was  made  for  its  simplicity.  Andrews' 
equation  of  state  exists  as  a deck  of  computer  cards,  and  ana- 
lytical expressions  describing  the  equilibrium  surfaces  are  not 
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readily  available.  However,  Andrews'  equation  of  state  does 
define  all  of  the  equilibrium  surfaces  of  alpha  and  epsilon  iron 
and  is  more  useful  than  the  present  equation  of  state  for  a 
variety  of  problems.  Except  for  the  assumptions  of  constant 
and  constant  ratio  of  Gruenisen  parameter  to  volume  for  the  epsi- 
lon phase,  the  present  equation  of  state  and  Andrews'  equation 
of  state  are  similar  in  principle.  In  practice,  the  two  ap- 
proaches differ  considerably.  For  example,  the  present  equation 
of  state  relies  on  one  set  of  isothermal  compression  data  to 
define  the  equilibrium  surfaces  while  Andrews  fitted  all  the 
available  thermodynamic  data  to  the  appropriate  derivatives  of 
a Helmholtz  potential  function  to  define  the  surfaces.  The  two 
equations  of  state  give  values  on  the  Hugoniot  for  volume  and 
temperature  at  130-kbar  stress  which  differ  by  less  than 
0.3  percent.  The  excellent  agreement  between  tne  equations  of 
state  was  expected. 

A.l.  Equilibrium  Surfaces 

The  specific  heat  at  constant  volume  is  assumed  constant 
following  the  approach  suggested  by  Johnson . This  approach 
is  sufficient  for  making  accurate  predictions  along  a Hugoniot. 
These  assumptions,  along  with  the  isothermal  compression  data  of 
alpha  and  epsilon  iron,  define  the  equilibrium  surfaces  for  the 
two  phases. 

The  primary  assumption  is  that  the  specific  heat  at  con- 
stant volume  is  independent  of  temperature  and  volume: 


constant 


(A.l) 
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This  assumption  restricts  E to  the  form, 

E = + f(V)  , 

where  £ (V)  is  an  arbitrary  function  of  specific  volume  V. 
Since  C„  is  constant.  Maxwell's  equations  require 


(A. 2) 


v/hich  gives 


P = r(V)T  + y{V)  = T - f (V)  . 


(A. 3) 


(A. 4) 


The  assumption  is  made  that 


nl  - c - - 

_ Jv  ■ SbEjy  ■ V 


= constant  , 


(A. 5) 


where  r is  t'\^  Iruneisen  parameter.  Then,  from  Eqs.  (A. 4)  and 


(A. 5) , we  find  that 


Cy  ^ = b Cy  = constant 


(A. 6) 


Now  any  path  on  the  thermodynamic  equilibrium  surface  can  be  used 
to  determine  y(V)  in  Eq.  (A. 4).  Therefore,  choosing  the  com- 
pression path  along  an  isotherm  is  a logical  choice  because  iso- 
thermal data  are  known  in  the  form  of  the  Murnaghan  equation. 

Each  phase  is  represented  by 


= fioi  ^ 


(A. 7) 
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L 


where  the  temperature  is  held  constant,  the  subscript  i 
indicates  the  particular  phase,  and  the  parameters  and  n^^ 
are  fitted  constants.  Combining  Eqs.  (A. 4),  (A.6) , and  (A. 7) 
results  in 


'V. 

1 


rv.  1 

-n . 

(T-Tq)  + Bq. 

1 

<1 

0 |h 

1  

^ - 1 

(A. 8) 


which  defines  any  P,V,T  path  on  the  equilibrium  surfaces. 

To  complete  the  definition  of  the  surfaces,  the  function 
f(V)  of  Eq.  (A. 2)  has  to  be  evaluated.  This  can  be  done  in  a 
consistent  way  with  the  use  of  the  following  identity: 


G*E-TS  + PV. 


The  use  of  subscript  i is  dropped  for  convenience. 
Eq.  (A. 9)  an  expression  for  S is  required.  Si^ce 


(A. 9) 

To  complete 


rasl  = h n 


(A. 10) 


and 


,/T 


(A. 11) 


Then 


S = Sn  - b C^(Vq-V)  + C„ln(T/T,,)  . 


(A. 12) 
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Using  the  identity. 


V =r^i 

L3pJt  » 


(A. 13) 


and  combining  Eqs.  (A. 2),  (A. 8),  (A. 9),  (A. 12),  and  (A. 13)  re- 


sults in 


- ■'0  >3  S * V„  - 1 . 


(A. 14) 


Integrating  eq.  (A. 14)  for  f (V)  and  substituting  in  Eq.  (A. 2) 


k ] 1 ''0"' 


results  in  the  expression. 


E = Eg  + C^(T-TqJ  - Tg  b C^(Vq-V)  + Vg  Bg 


(A. 15) 

Equations  (A. 8),  (A. 12),  and  (A. 15)  define  the  phase  1 

and  phase  2 equilibrium  surfaces  within  initial  values  Vg^, 

V02,  Egj^,  Eg2,  Sgj^,  snd  Sq2.  The  value  for  Vgj^  for  iron 

is  easily  measured  while  Vg2  has  been  inferred  from  x-ray  meas- 
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urements  of  coexisting  phases.  The  value  for  Sg^^  can  be 
arbitrarily  chosen  while  Sg2  is  fixed  by  the  difference  in 
entropy  between  phases  which  can  be  inferred  from  the  measured 
slope  of  the  phase  line  or  obtained  from  a thermodynamically 
consistent  empirical  equation  of  state  which  is  fitted  to  all 
types  of  experimental  data.  ' The  magnitude  of  Egj^  can  be 
arbitrarily  chosen  while  Eg2  must  be  specified  by  some  addi- 
tional relationship.  A relationship  which  makes  the  connection 
between  the  equilibrium  surfaces  is: 
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G2i(P*,T*)  = 0 , (A. 16) 


where  62^  is  the  difference  between  Gibbs  energies  of  the 

phases,  (P*,T*)  is  the  equilibrium  transformation  point.  This 

equilibrium  condition  with  Eq.  (A. 9)  determines  the  value  of 

Eq2*  Use  of  Eq.  (A. 16)  to  determine  the  equilibrium  surface 

provides  no  assurance  that  the  state  point  in  any  real  process 

will  actually  lie  on  such  a surface  (cf  Chapter  5). 

The  values  for  E^j^,  and  8^2  st  i atmosphere 

pressure  and  a temperature  of  295°K  were  chosen  to  be  the 

27  26 

same  as  those  used  by  Andrews.  ' The  value  for  Eq2  was 
determined  from  the  relation  G22^(116  kbar,332°K)  = 0.  The 
initial  value  of  C^2  classical  value  of  3R  where  R 

is  the  gas  constant.  The  initial  values  used  to  define  the 
equilibrium  surfaces  are  listed  in  Table  A.l. 


TABLE  A.l.  — initial  values 
(T=295  K,P=0) 


Variable 

Alpha  phase 

Epsilon  phase 

Vq 

1 

0.1272 

0.1208  cm^/gm 

7.932  X 10"^ 

1.667  X 10”^  Mbar  cm^/gm 

Sq 

4.863  X 10"® 

5.323  X 10~®  Mbar  cm^/gm  deg 

n 

5.917 

5.102 

»0 

0.275 

0.325  Mbar 

r 

1.735 

2.0 

S 

4.447  X 10“® 

4.466  X 10"®  Mbar  cm^/gm  deg 
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A. 2.  Approximations  of 
Potential  Functions 

It  is  often  found  convenient  to  use  simple  approximations 
for  the  differences  between  the  internal  energies  of  the  two 
phases  and  also  the  difference  between  the  Gibbs  energies.  To 
demonstrate  their  accuracy,  the  approximations  are  compared  to 
the  more  accurate  equation  of  state  results.  It  is  important  to 
keep  in  mind  that  these  differences  are  between  values  that  lie 
on  the  equilibrium  surfaces.  The  first  useful  approximation  is 
for  G2  - Gj^.  Since 

d(G2-G^)  = (V2-V^)dP  - (S2-Sj^)dT  ; (A. 17) 

if  V2  - and  S2  - are  constant,  then 

G2  - Gi  * (V2-Vj^)  (P-P*)  - (S2-Sj^)  (T-T*)  , (A.  18) 

the  difference  between  internal  energy  of  the  two  phases,  is 

E2  - = T*(S2-Sj^)  - P*(V2-V^)  . (A.  19) 

For  the  approximate  calculations,  the  following  values 

3 -7  3 

are  used:  V2-Vj^  = -0.00604  cm  /gm,  ^2-5^  ~ 3.67  x 10  Mbar  cm  / 

gm,  and  P*  = 0.116  Mbar,  T*  = 332®K.  Figure  A.l  shows  the 
result  from  the  approximate  relation  for  the  Gibbs  energy  differ- 
ences and  those  from  the  more  exact  equation  of  state  for  the 
mixed  phase  Hugoniot  of  iron.  The  agreement  between  Eq.  (A. 18) 
and  the  equation  of  state  results  is  excellent.  The  approximate 


f 


‘I 

\ 
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value,  E2-Ej^  = 8.22  x 10~  Mbar  cm^/gm,  agrees  within  2 per- 
cent of  the  equation  of  state  results  over  the  stress  range  of 
130  to  200  kbar. 
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APPENDIX  B 

STREAK  CAMERA  CALIBRATION 

The  Cordin  132  streak  camera  writing  speed  was  calibrated 
by  recording  the  output  of  a crystal-controlled  oscillator  dis- 
played on  the  face  of  an  oscilloscope.®^  The  dynamic  spatial 
resolution  for  the  camera  was  measured  by  recording  the  image  of 
a Ronchi  ruling  for  different  camera  writing  speeds.  The  details 
of  the  calibration  experiments  are  given  below. 

B.l.  Writing  Speed 

A schematic  of  the  calibration  setup  is  given  in  Fig.  B.l. 
The  equipment  used  for  the  calibration  experiments  consisted  of 
the  streak  camera,  the  synchronizer  circuit,  a Hewlett-Packard 
5512A  electronic  counter,  a resistance-capacitor  (RC)  pulse 
power  supply,  a Tektronix  535  oscilloscope,  and  a General  Radio 
1213A  oscillator.  The  oscillator  was  stable  to  within  +4  hz 

7 

out  of  10  hz  for  24  hours.  The  frequency  scale  of  the 
oscillator  was  accurate  to  +3.4  hz  out  of  10®  hz.  The  elec- 
tronic counter  was  accurate  to  one  count. 

The  electronic  equipment  was  connected  as  shown  in 
Fig.  B.l  and  the  camera  was  focused  on  the  face  of  the  oscillo- 
scope's cathode  ray  tube.  The  camera's  mirror  was  rotated  at  a 
constant  speed  and  then  a switch  on  the  camera  controls  was 
closed  which  opened  the  camera  shutter  and  sent  a pulse  to  the 
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RC  POWER 
SUPPLY 


CAMERA 

SYNCHONIZER 

CIRCUIT 


FIGURE  B.1  SCHEMATIC  OF  CALIBRATION  TECHNIQUE  FOR  THE  STREAK  CAMERA. 
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synchronizer  circuit.  The  synchronizer  circuit  assures  that  the 
rotating  mirror  is  in  proper  position  to  record  on  the  film  and 
supplies  a pulse  to  trigger  the  RC  power  supply.  The  power  sup- 
ply applied  -80V  pulse  to  the  cathode  ray  tube  which  brightens 
the  scope  trace  for  a few  milliseconds.  A 100  Khz  signal  from 
the  oscillator  was  continuously  applied  to  the  vertical  input  of 
the  oscilloscope  without  sweeping  the  scope  in  time.  A vertical 
oscillating  line  image  of  the  oscillator's  signal  was  displayed 
on  the  scope  face  and  then  recorded  dynamically  by  the  streak 
camera.  A sinusoidal  wave  results  on  the  streak  camera  film  with 
accurately  known  peak-to-peak  time  intervals. 

The  camera  writing  speed  is  given  by 

W = F • R (B.l) 

where  W is  the  writing  speed,  F is  the  calibration  parameter, 
and  R is  the  rotation  speed  of  the  camera's  mirror.  The  re- 
sults of  the  calibrations  are  given  in  Table  B.l. 

TABLE  B.l. — Calibration  experimental  results 


Experiment 

Number 

R 

{ rev/sec) 

W 

(mm/psec) 

F 10"^ 

(mm) 

Cal-1 

250 

1.0075  + .0025 

4.030  + .001 

Cal-2 

250 

1.0075  + .0025 

4.030  + .001 

Cal-3 

1000 

4.04  + .02 

4.040  + .02 
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The  record  for  experiment  Cal-3  was  faint  and  could  not 

be  read  precisely  as  the  given  error  indicates,  but  the  results 

are  in  good  agreement  with  Cal-1  and  Cal-2.  The  calibration 

parameter  F should  be  constant  for  revolution  speeds  up  to 

62 

5000  hz  since  beryllium  mirrors  do  not  physically  distort  for 
these  rotation  speeds. 

B.2.  Dynamic  Spatial  Resolution 

The  practical  definition  of  the  spatial  resolution  was 
taken  as  the  number  of  lines/mm  distinguishable  on  the  film  in 
a simulated  experimental  setup.  The  dynamic  spatial  resolution 
of  the  camera  was  measured  by  dynamically  recording  the  image  of 
a Ronchi  ruling.  The  glass  ruling  plate  was  located  in  the  test 
chamber  so  the  camera  would  view  it  in  an  experimental  situation. 
The  camera  was  focused  on  the  ruling  with  the  projection  of  the 
slit  plane  perpendicular  to  the  rulings.  The  light  for  the 
camera  was  provided  by  an  exploding  wire  light  source  located 
just  behind  the  ruling.  The  camera  was  run  at  speeds  up  to 
15  mm/usec.  The  dynamic  spatial  resolution  for  this  range  of 
writing  speeds  was  > 30  lines/mm. 
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APPENDIX  C 


ERROR  ANALYSIS 


This  appendix  outlines  the  error  analysis  used  for  this 
work.  The  qeneral  treatment  of  errors  is  given  in  Section  C.l, 
and  a specific  example  of  the  error  in  the  measurements  of  free 
surface  velocity  is  treated  in  Section  C.2. 

C.l.  General  Treatment  of  Errors 

The  error  in  a quantity  is  the  uncertainty  in  its  meas- 
urement. The  value  of  a quantity  may  depend  on  two  or  more  other 
parameters.  Each  of  these  parameters  is  uncertain  to  some  degree. 
To  estimate  the  total  uncertainty  in  the  quantity  requires  a 
prescription  for  calculating  the  error  which  includes  the  propa- 
gation of  errors.  The  present  error  analysis  gives  such  a 
prescription  which  depends  on  obtaining  analytical  expressions 
of  the  error  in  terms  of  known  parameters.  The  analysis  is  based 
on  a book  by  Yard  ley  Beers. 

The  errors  are  reported  in  the  short  liand  form  "f  + i" 
where  r is  tlie  average  deviation  of  f.  The  total  magnitude 
of  the  error  is  just  2i . 

Errors  may  be  random,  as  in  the  case  when  a measurement 
is  repeated  several  times.  The  results  will  be  distributed 
around  a "most  probable"  value  which  is  assumed  here  to  be  the 
average  value.  For  random  errors  there  is  a possibility  of 
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compensation  among  the  various  contributions.  It  is  expected 
that  the  total  error  will  be  algebraically  less  than  the  sum  of 
the  separate  contributions.  A logical  way  of  adding  the  sepa- 
rate contributions  is  to  take  the  square  root  of  the  sum  of  their 
squares  (hereinafter  called  SRSS)  which  does  have  the  compensat- 
ing property.  The  rule  for  combining  random  errors  is 


where  f - F(x,y)  and  x,y  are  the  independent  measured  param- 
eters. The  function  f can  depend  on  any  number  of  parameters. 

The  choice  of  two  parameters  for  Eg.  (C.l)  was  made  only  for 
illustrative  purposes. 

Errors  may  be  systematic  where  all  the  individual  values 
are  in  error  by  the  same  amount.  These  errors  cannot  be  esti- 
mated by  repeated  measurements.  The  systematic  error  often  comes 
from  the  uncertainty  in  the  calibration  of  the  instrument  used 
to  measure  the  parameter.  These  errors  combine  algebraically 
because  compensation  among  the  various  contributions  is  not 
likely.  Therefore, 

^f  " 3x  ^x  ^ ^y  * 

When  both  kinds  of  errors  are  encountered  they  are  com- 
bined by  taking  the  SRSS  of  both  separate  contributions  to 
give  a total  estimate  of  error  for  that  particular  parameter. 
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From  a set  of  measurements  of  a parameter  the  best  value 


available  is  the  average  value.  The  error  associated  with  this 
average  value  ia  not  the  average  of  the  deviations  from  this 
average  value  but  rather  the  average  of  these  deviations  divided 
by  the  square  root  of  the  number  of  measurements.  In  other 
words,  the  precision  improves  in  proportion  to  the  square  root 
of  the  reciprocal  of  the  number  of  measurements  in  the  sample 
set. 


C.2.  Error  in  Free  Surface 
Velocity  Measurements 

The  experimental  technique  used  to  measure  the  free  sur- 
face motion  was  to  photographically  record  with  a streak  camera 
the  changes  in  the  intensity  of  reflected  light  from  the  tilted 
mirrors  which  were  placed  on  the  surfaces  of  the  test  samples. 

For  this  experimental  technique  the  free  surface  velocity  is  a 
function  of  five  different  parameters.  (See  Chapter  2.)  Each 
parameter  is  unknown  by  a certain  amount.  Each  individual  uncer- 
tainty in  the  five  parameters  contributes  to  the  overall  uncer- 
tainty in  the  free  surface  velocity.  Determination  of  the  total 
uncertainty  is  the  goal  of  this  appendix. 

The  free  surface  velocity  due  to  the  plastic  I shock  is 


given  by 


tan  g 

tan  Y - tan  o) 


(C.3) 


where  u^^  is  the  free  surface  velocity  due  to  the  plastic  I 
shock,  u is  the  camera  speed,  a is  the  tilted  mirror  angle. 
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Mj  is  the  magnification  (image/object),  y is  the  angle  of  the 
trace,  and  w is  the  wave  tilt  angle  on  the  film.  The  total 


The  individual  errors  for  camera  speed,  magnification, 
mirror  angle,  trace  angle,  and  wave  tilt  angle  contain  both  ran- 
dom and  systematic  errors.  The  random  errors  were  estimated  by 
taking  the  deviations  from  the  average  of  a number  of  measure- 
ments of  each  parameter.  The  systematic  error  was  usually  taken 
as  equal  to  one-half  of  the  smallest  division  of  the  instrument 
used  to  measure  the  parameter.  The  total  error  for  each  param- 
eter was  then  found  by  taking  the  SRSS  of  the  systematic  and 
random  errors. 


C.2.1.  Sources  of  Various  Errors 

The  performance  of  accurate  experiments  requires  knowl- 
edge of  the  source  and  magnitude  of  the  various  errors.  The 
purpose  of  this  subsection  is  to  identify  the  sources  of  error 
in  the  free  surface  velocity  measurement  and  to  determine  which 
are  the  important  errors. 

The  error  in  the  camera  writing  speed  was  determined 
from  the  error  associated  with  the  calibration  constant  (see 
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Appendix  B)  and  the  uncertainty  of  +1  count  in  the  counter 
reading  which  determined  the  camera’s  rotor  speed. 

The  total  error  in  the  magnification  was  determined  by 
the  length  of  the  precision  grid  used  and  the  uncertainty  of 
+0.05  mm  in  measuring  the  image  of  the  grid  with  a precision 
eyepiece. 

The  total  error  in  the  mirror  angle  a was  determined 
by  the  scatter  in  the  readings  and  the  uncertainty  of  +1  min  of 
arc  in  measuring  the  angle  with  a Shadowgraph  instrument. 

The  error  in  determining  the  trace  angle  y and  the  tilt 
angle  tii  is  more  difficult  to  ascertain.  They  require  an  under- 
standing of  the  camera  operation,  experimental  design,  and 
record  reading  technique.  A detailed  treatment  of  this  error  is 
given  in  Section  C.2.2. 

Table  C.l  summarizes  the  total  contributions  to  the 
error  in  u^^  for  a typical  experiment.  The  largest  and  most 
important  error  is  in  the  trace  angle  y. 

TABLE  C.l. — Various  contributions  to  the  error  in  u^^ 


% 

S 

"u) 

^c 

«f 

sin  2ot 

2 

cos  y{tan  y - tan  co) 

5 

cos  u)(tan  y - tan  oj) 

0.006 

0.005 

0.008 

0.028 

0.003 
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C.2.2.  Errors  in  Determining 
the  trace  Angle  and  the 
Wave  Tilt  Angle 

It  was  shown  above  that  the  error  in  determining  the 
trace  angle  on  the  film  is  the  largest  contributor  to  the  total 
error  in  Therefore,  a detailed  discussion  on  the  syste- 

matic part  of  this  error  is  given  here.  The  estimates  of  the 

systematic  errors  in  determining  the  trace  angles  on  the  camera 

64 

records  are  similar  to  those  of  Dubovik.  There  are  four  sepa- 
rate contributions  to  this  error:  (1)  the  film  alignment  in  the 
reader  instrument,  (2)  the  fuzziness  of  the  trace,  (3)  the  reader 
instrument  accuracy,  and  (4)  the  spatial  resolution.  Each  is 
treated  as  independent  of  the  other,  which  is  not  strictly  true. 

It  is  a reasonable  assumption  to  make  since  the  sole  purpose  is 
to  obtain  a reasonable  estimate  of  the  error. 

The  first  contribution  to  the  systematic  error  in  y 
is  due  to  the  alignment  of  the  film  parallel  to  the  time  axis 
for  reading.  This  error  is  <0.001  rad  based  on  a graphic 
accuracy  of  +0.1  mm  per  100  mm  which  can  be  realized  in 
practice . 

The  second  and  largest  contribution  to  the  error  in  y 
is  due  to  the  fading  and  fuzziness  of  the  trace,  which  depends 
on  the  camera  slit  width  and  response  of  the  film  to  different 
exposure  times.  The  recording  of  the  motion  of  an  opaque  half- 
plane moving  at  a constant  speed  parallel  to  the  slit  plane 

(along  the  length  of  the  slit)  results  in  a cut-off  of  the  ligh*-  j 

1 

to  the  camera.  The  cut-off  is  not  recorded  as  a sharp  loss  of  ; 

light  because  of  the  finite  width  of  the  slit.  Rather,  it  is  a ^ 

1 

)i 
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line  of  some  width  made  up  of  many  thin  lines  each  at  different 
optical  densities  and  at  an  angle  y to  the  space  axis  of  the 
film.  The  optical  density  of  the  broad  line  varies  from  high  to 
low  density  in  the  direction  of  the  camera  sweep  (e.i.,  increas- 
ing time) . The  geometry  of  the  envelope  of  lines  is  shown  in 
Fig.  C.l. 

An  analytical  expression  for  the  error  can  easily  be 
derived  from  the  geometry,  resulting  in 


e = AB  sin  y/tS-AB  cos  y) 

where  AB  is  the  camera  slit  width,  y is  the  angle  of  the 


(C.5) 


trace,  and  S is  the  trace  length.  The  typical  systematic  error 
is  0.018  rad  for  AB  = 0.05  mm,  y = 0.786  rad,  and  S = 2 mm. 

The  third  contribution  to  the  error  is  that  due  to  the 
accuracy  of  the  Vanguard  analyzer  in  determining  y.  By  assuming 
the  error  can  be  assigned  to  one  reader  variable,  then  Eq.  (C.5) 
is  applicable.  The  slit  width  AB  in  Eq.  (C.5)  is  replaced  by 
the  smallest  distance  measurable  by  the  Vanguard  analyzer 
(0.005  mm).  A typical  error  is  0.002  rad  for  y = 0.786  rad 
and  S » 2 mm. 


The  fourth  contribution  to  the  error  is  due  to  the  spa- 
tial resolution.  The  smallest  distinguishable  spatial  distance 
on  the  film  was  about  0.03  mm  for  the  camera  in  a simulated 
experimental  setup.  (See  Appendix  B.)  Treating  this  uncertainty 
as  independent  of  the  others,  an  expression  for  this  error 


becomes 


0.03  mm 


• cos  y . 


(C.6) 
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The  typical  contribution  for  this  error  is  0.011  rad  when 
S = 2 nun  and  y = 0.786  rad. 

The  typical  errors  were  given  above  for  each  separate 
case  to  illustrate  which  was  the  largest  and  dominating  error. 

It  is  obvious  by  comparing  the  typical  contributions  given  for 
each  factor  above  that  the  largest  errors  are  due  to  the  fuzzi- 
ness of  the  trace  and  to  the  spatial  resolution  of  the  camera. 

It  should  be  emphasized  that  the  above  analysis  for  the 
error  in  y is  a pessimistic  one.  The  analysis  was  limited  to 
the  error  between  two  data  points  while  the  actual  measurement 
of  the  angle  typically  resulted  in  15  data  points  which  were 
least-squares  fitted  to  a straight  line.  The  least-squares 
fitting  of  the  data  reduces  the  measuring  error.  In  addition, 
placing  the  Vanguard  analyzer's  hairlines  on  film  coordinates  of 
about  the  same  optical  density  reduces  the  trace  width  to  less 
than  the  slit  width.  To  compensate  for  these  factors,  personal 
judgment  must  be  used.  A reduction  of  one-half  of  the  error  due 
to  the  slit  width  error  in  Eq.  (C.5)  was  used. 

The  same  methods  used  in  estimating  errors  in  the  trace 
angle  apply  and  were  used  for  the  errors  in  the  wave  tilt  angle 
u>. 

il 
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APPENDIX  D 

PERMANENT  REGIME  WAVE  SOLUTION 

57 

A nonsteady  shock  will  approach  a steady  shock  in  time. 

A solution  for  the  steady  shock  is  given  here  based  on  the  const i- 

20 

tutive  relation  of  the  Horie-Duvall  model. 


D.l.  General  Solution 

The  differential  flow  equations  in  Eulerian  coordinates 

are: 


(D.l) 


« du  3P 
P dt  ■ “ 3x  ' 


(0.2) 


dE 

dt 


P s 


1 

V • 


(0.3) 


The  total  time  derivative  is  (d/dt)  = (3/3t)  + u(3/3x).  The 
steady  state  condition  is  that  3/3t  = 0 which,  applied  to  the 
flow  equations,  results  in 


pu  = m = constant  , (D.4) 

P + mu  = constant  , (D.5) 
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(D.5)  results  in  a useful  equation  which  defines  the  compression 
path  as  a straight  line  in  the  P-V  plane  connecting  the  initial 
and  final  P-V  states.  The  equation  is 

P - P*  = -ro^(V-V*)  . (D.7) 

Another  useful  relation  is  obtained  from  Eqs.  (D.l)  and  (D.2); 


dV  1 dP  . 

cHc  ~ Hx  * ' 

m 

In  the  mixed  phase  region  the  extensive  parameters  V 
and  E are  defined  as  functions  of  P,  T,  and  f such  that; 

V(P,T,f)  « V^(P,T)  + fCV2(P,T)  - Vj^(P,T)]  , ( 

E(P,T,f)  = Ej^(P,T)  + f[E2(P,T)  - Ej^(P,T)]  . I 

Since  P,  T,  and  f are  implicit  functions  of  x,  Eq.  (D.9) 
and  (0.10)  result  in 


(0.8) 


[avl 

dP  , 

'avl 

dx 

[apj 

aTj 

•is  1 

'3E 

dP 

^3e1 

“ 1 

P. 

T,f  ^ 

[jtJ 

r avl  dT  . f wl  df 

LHp.f  3* 


(D.9) 

(0.10) 


(D.ll) 


(0.12) 


Using  Eqs.  (0.8),  (0.9),  and  (D.IO)  and  the  thermodynamic  iden- 
tities for  specific  heat  Cp,  compressibility  K^,  and  thermal 
expansion  8,  then  Eqs.  (0.11)  and  (0.12)  become: 
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dP 

33?  ~ 


VBf  dT  , V \ 

^ 33?  33? 


T,f 


m 


9 


dT 

33?  = 


dp  r _p 


+ TVSf 


=P,f  - '’''Sf 


(D.13) 


(D.14) 


Equations  (D.13)  and  (D.14)  require  only  an  equation  of  state  and 
a relation  for  df/dx  to  solve  for  the  stress  and  temperature 
distribution  in  a steady  wave. 

For  iron  the  second  term  on  the  right  hand  side  of 
Eq.  (D.13)  is  about  60  times  greater  than  the  first  term. 

Steady  profile  calculations  with  this  first  term  neglected  (for 
example  see  Fig.  (6.2)  are  called  temperature  independent. 


D.2.  Rate  Equation 

20 

Hone  and  Duvall  assumed  a relation  for  df/dt  (see 
Chapter  6)  that  for  steady  waves  can  be  written  as 


df  _ f®'^  - f 
33? T~ 


(D.15) 


27  29 

Andrew  later  showed  that  by  defining  t as  the 

constant  energy  and  constant  volume  equlibration  time,  the  rate 
equation  can  be  directly  related  to  the  difference  in  Gibbs 
energies.  For  steady  waves  this  rate  equation  becomes: 


U 


df 

33?  ■ " 


J G 


21 


IA|  T. 


(D.16) 


) 

1 
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For  small  variations  from  equilibrium  the  determinant 


1 

SfV 

AV 

A1  = 

PVX,p^f-TVef 

Cp^f-PV3f 

AE 

AV 

-AS 

0 

= Aa(eVA-7-VK^AS)  + AVAS(TV3-PVK^)  + ( AV)  ^ (PV3-Cp) 

(D.17) 

where  AS  = S2-Sj^  and  AV  = 

The  transformation  Jacobian  J was  previously  evaluated 
by  Hayes^^  as 


* j[m!]  * -‘^P,£™S,(  • '■’•18) 

Equation  (D.16)  provides  the  necessary  relation  defining 
df/dx.  Equations  (D.13),  (D.14),  {D.17),  and  (D.18),  and  an 
equation  of  state  give  the  stress  and  temperature  distribution 
in  the  steady  shock. 
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